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Dear  Mr.  Smith, 

Attached  please  find  the  FY98  Annual  Interim  Technical  Report  for  the  SERDP  Project  1094, 
“Environmental  Impacts  to  the  Chemical  Signature  Emanating  from  Buried  UXO.”  The  authors  hope 
that  this  documentation  will  provide  sufficient  technical  documentation  of  progress  for  the  first  year  of 
this  exciting  project.  We  feel  that  the  findings  will  provide  the  needed  information  to  scientific  and 
operational  groups  to  further  the  development  of  chemical  sensing  for  detection  of  UXO.  If  there  are 
any  questions,  please  contact  myself  at  your  convenience. 


Sincerely, 


Exceptional  Service  in  the  National  Interest 


1 .0  Project  Background 

The  goal  of  locating  buried  UXO  and  landmines  is  a  significant  challenge  to  science  and 
technology.  Technology  development  efforts  are  under  way  to  use  chemical  sensors  that 
can  discriminate  inert  ordnance  and  clutter  from  live  munitions  that  continue  to  be  a 
threat.  However,  the  chemical  signature  is  affected  by  multiple  environmental 
phenomena  that  can  enhance  or  reduce  its  presence  and  transport  behavior,  and  can  affect 
the  distribution  of  the  chemical  signature  in  the  environment.  For  example,  the  chemical 
can  be  present  in  the  vapor,  aqueous,  and  solid  phases.  The  distribution  of  the  chemical 
among  these  phases,  including  the  spatial  distribution,  is  key  in  designing  appropriate 
detectors,  e.g.,  gas,  aqueous  or  solid  phase  sampling  instruments.  A  fundamental 
understanding  of  the  environmental  conditions  that  affect  the  chemical  signature  is 
needed  to  describe  the  favorable  and  unfavorable  conditions  of  a  chemical  detector  based 
survey  to  minimize  the  consequences  of  a  false  negative.  Figures  1  and  2  show  the 
principal  subsurface  phenomena  and  the  principal  surface  boundary  conditions  that  are 
being  explored  in  this  project. 


Chemical  Detection 
System 


Ground  Surface^ 


Figure  1 .  Subsurface  Environmental  Fate  and  Transport  Phenomena 


2.0  Objective 

The  objective  of  this  project  is  to  develop  a  validated  subsurface  transport  model  that  can 
be  used  to  predict  the  spatial,  temporal,  and  phase  specific  concentration  of  chemical 
signature  molecules  derived  from  shallow  UXO  and  buried  landmines  under  the 
influence  of  specific  environmental  conditions.  Other  government  programs  are 
developing  chemical  detector  platforms  that  can  provide  a  separate  unique  signal  to 
classify  subsurface  objects  identified  with  existing  geophysical  systems.  It  is  estimated 
that  eleven  million  acres  of  land  needs  assessment  to  identify  subsurface  UXO  -  with 
costs  estimated  to  be  about  $1 .4M/acre.  The  ranges  where  UXO  can  be  found  are 
distributed  throughout  the  country  where  environmental  conditions  vary  significantly.  It 
is  the  hypothesis  of  this  project  that  these  environmental  conditions  will  have  a 
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significant  impact  on  the  transport  of  chemical  signature  molecules  from  subsurface 
UXO  and  buried  landmines  to  the  surface  before  presentation  to  a  chemical  detector 
system. 

If  through  this  systems  analysis,  one  can  show  the  ranges  and/or  combinations  of 
environmental  parameters  that  improve  the  transport  of  chemical  signature  molecules  to 
the  chemical  detector  system,  and  conversely,  those  that  constrain  this  movement,  end- 
users  seeking  to  locate  buried  UXO  will  be  better  positioned  to  understand  the  merits  and 
limitations  when  looking  to  deploy  the  chemical  detector  technology. 


Figure  2.  Soil  Surface  Boundary  Conditions 


3.0  Technical  Approach 

This  project  is  a  cooperative  effort  involving  Sandia  National  Laboratories,  US  Army 
Corps  of  Engineers  -  Cold  Regions  Research  and  Engineering  Laboratory,  the  University 
of  California  Riverside  and  Lawrence  Berkeley  Laboratory.  This  project  has  been  co¬ 
funded  by  the  DARPA  program  for  the  chemical/biological  detection  of  explosives  in 
landmines.  The  project  is  divided  into  five  major  tasks. 

Task  1:  Model  Development/Utilization  -  The  task  involves  completion  of  a  sensitivity 
analysis  of  known  input  parameters  in  a  one-dimensional  analytical  contaminant  transport 
model,  expanding  this  model  to  assess  two-dimensions  to  explore  the  surface  area 
footprint  from  buried  UXO,  and  modifying  an  existing  numerical  simulation  code 
(TOUGH2)  (surface  boundary  conditions,  multiple  chemical  components,  liquid 
diffusion)  for  use  as  the  complete  systems  analysis  tool.  Inverse  modeling  will  be  used  to 
assess  input  parameter  sensitivity  and  as  a  tool  for  the  design  of  laboratory  validation 
experiments  in  task  three. 

Task  2:  Fundamental  Properties  -  This  task  involves  the  measurement  of  specific 
transport  parameters  currently  not  available  in  the  literature  for  explosive  signature 
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molecules.  These  include  temperature  dependent  water  solubility,  vapor-solid  sorption  as 
a  function  of  soil  moisture  content  and  source-term  emission  rates. 

Task  3:  Laboratory-Scale  Experiments  -  This  task  involves  a  laboratory  validation  study 
that  will  confirm  the  most  critical  parameters  included  in  the  simulation  model  through 
soil  column  transport  experiments. 

Task  4:  Operational  Strategy  -  This  task  will  utilize  the  simulation  tools  to  assess  the 
impacts  of  environmental  conditions  on  the  transport  of  chemical  signature  molecules 
from  shallow  UXO  and  support  end-user  queries  on  the  utility  of  chemical  sensor 
platforms  for  the  classification  stage  in  the  identification  of  true  unexploded  ordnance. 

Task  5:  Ordnance  Source  Term  -  This  task  was  added  in  August  1998  to  evaluate  the 
source  term  from  actual  ordnance  items. 

4.0  Project  Accomplishments 

Task  One  -  Model  Development/Utilization 

The  environmental  fate  and  transport  of  organic  chemicals  including  volatilization  and 
leaching  losses  has  been  used  to  explore  the  distribution  of  agricultural  pesticides  in  soils 
(Mayer  et  al.  1974,  Farmer  et  al.  1980,  and  Jury  et  al.  1980).  These  models  were 
primarily  intended  to  simulate  specific  circumstances.  However,  Jury  et  al.  (1983,  1984a, 

1984b,  1984c)  developed  and  validated  a  general  screening  model  (Behavior  Assessment 
Model,  BAM)  that  included  volatilization,  leaching,  and  degradation  to  explore  the  major 
loss  pathways  of  agricultural  pesticides  as  a  function  of  specific  environmental 
conditions.  The  Behavior  Assessment  Model  was  adapted  for  evaluation  of  chemicals  in 
buried  soils  and  has  been  termed  the  Buried  Chemical  Model  (BCM)(Jury  et  al.,  1990). 

The  model  simulations  can  be  used  to  assess  the  behavior  of  different  chemicals  under 
particular  environmental  conditions,  but  are  not  intended  to  predict  a  definitive 
concentration  distribution  in  the  field.  As  such,  the  predictions  from  the  screening  model 
are  only  an  indication  of  expected  conditions. 

A.  Phase  Partitioning  Phenomena 

The  formulations  of  the  BAM  and  BCM  models  begin  by  defining  phase  partitioning 
phenomena.  These  are  valuable  in  that  they  can  express  the  total  concentration  of  a 
chemical  in  the  gas,  aqueous  and  sorbed  phases.  The  total  concentration  is  expressed  as 

Cj.  =  PhCg  +  dCi  +ciCq  [1] 

where  Cs  is  the  concentration  sorbed  to  the  soil.  Cl  is  the  solute  concentration  in  the 
aqueous  phase,  and  Cg  is  the  gas  phase  concentration.  In  addition.  Jury  (1983)  shows 
how  equation  [1]  can  be  rewritten  in  terms  of  one  of  the  variables  alone 

C,  =  R,C,  =  R,C,  =  [2J 

where 

^.v  =  A +  ^ 
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R,  =  +  ^  +  ,  and 

[4] 

Kj  e 

[5] 

are  the  solid,  liquid  and  gas  phase  partition  coefficients,  respectively. 

In  their  evaluation  of  vapor  phase  transport  in  soils,  Ong  et  al.,  1992,  added  vapor-solid 
sorption  such  that  equations  [1]  and  [2]  become 

Cj-  —OCi  +ciCq  +C^KjPi,  -hCQKgQPi, 

[6] 

and 

Ct  =  Rs^S  ~  ~  ^G^G  ~  ^SG^SG 

[7] 

where 

Rj  =  Pi,Kj  +  6  +  +  Pi^KfjKgQ 

[8] 

n  K,  0 

^G=Ph  +  „  +^  +  PhKsG 

[9] 

P  r,  \  ^  \  1  Pi>^H^SG  j 

RsL=Ph+  r.  ^ 

^d  ^d 

[10] 

p^Kj  6  a 

—  +  H-  + 

IT  IT  ^ " 

^H^SG  ^H^SG  ^SG 

[11] 

are  the  liquid,  gas  and  solid-liquid  and  solid-gas  phase  partition  coefficients,  respectively. 
This  formulation  introduces  a  new  term,  Ksg  that  is  a  function  of  the  overall  vapor 
partition  coefficient  (K'd),  which  is  highly  dependent  on  the  soil  moisture  content.  Ksg  is 
defined  as  (Ong  et  al.,  1992) 


Below  four  monomolecular 
function  described  by 


w 


layers  of  water  coverage  on  soils,  K'j  is  an  exponential 


A  =  \og{K\ ) 

A  =  {A,-  +  p{w) 


[12] 


[13] 

[14] 
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[15] 


a  = 
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^4-/g 
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J3(w)  =  log  • 
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Above  four  monomolecular  layers  of  water  coverage  on  soil,  K'd  is  a  function  of  Kh  and 
is  described  by 


[16] 


f 

A  =  log 

V 


El. 

Kh 


+ 


w  '] 


KhYPw  J 


Ong  et  al,  1990,  characterize  the  vapor-solid  partitioning  in  this  region  as  being 
controlled  by  Henry’s  Law  Constant  (HLC).  This  is  because  the  vapor  must  first 
partition  into  the  soil  water  prior  to  partitioning  onto  the  soil  particle.  The  moisture 
content  at  four  monomolecular  layers  of  soil  water  is  a  function  of  the  soil  specific 
surface  area  and  is  described  by 


[17] 


W4 


MA^A„  ^ 


The  specific  surface  area  (Sa)  of  soils  range  from  10  mVg  for  sand  to  lOO’s  mVg  for 
some  types  of  clay.  Figure  3  shows  the  correlation  of  soil  specific  surface  area  to  the 
moisture  content  at  four  monomolecular  layers  of  water.  The  proportion  of  clay  in  soils 
strongly  influences  the  soil  specific  surface  area.  This  is  due  to  the  small  size  of  the  clay 
soil  particles. 


[18] 


Figure  3.  Soil  Water  Content  at  4  Monomolecular  Layers  as  a  Function  of  Soil  Specific  Area 
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In  their  evaluation  of  toluene  and  trichloroethene,  K'j  increased  about  10'*  from  the  point 
of  four  monomolecular  layers  to  oven  dry  soil  moisture  contents.  The  impact  of  the 
vapor  solid  partitioning  is  significant  at  moderately  low  soil  moisture  contents.  Figure  4 
shows  how  the  relative  soil  gas  concentration  can  decline  rapidly  as  the  soil  moisture 
content  declines.  Note  that  the  soil  moisture  contents  at  four  monomolecular  layers  are  at 
levels  that  are  not  unusually  low,  and  are  typically  observed  in  soils  after  precipitation 
and  drainage  events  have  occurred.  With  an  extended  absence  of  precipitation,  the 
surface  soil  moisture  approaches  the  extremely  dry  region.  It  is  a  potential  that  this  dry 
layer  could  act  as  a  barrier  to  vapor  emission  from  soils  and  be  a  preconcentrator  of 
analyte  signal.  Whether  this  process  is  fully  reversible  with  the  addition  of  water  is 
unknown,  and  may  represent  an  opportunity  for  enhancement  of  the  chemical  signature. 


0  0.05  0.1  0.15  0.2  0.25 

Soil  Moisture  Content  (g/g) 

Figure  4.  Relative  Soil  Gas  Concentrations  Impacted  by  Low  Soil  Moisture  Content 

In  order  to  understand  this  phenomenon  more,  a  Monte  Carlo  simulation  approach  was 
employed.  This  method  specifies  a  statistical  distribution  function  for  various  input 
parameters,  then  randomly  selects  particular  values  for  each  input  variable  within  the 
specified  statitical  distribution,  and  tracks  the  output  of  one  or  more  variables.  The 
following  input  assumptions  and  distributions  were  defined  for  the  soil  bulk  density,  soil 
particle  density,  soil  moisture  content,  soil  water  partition  coefficient,  soil  specific 
surface  area  and  soil  vapor  partition  coefficient  at  oven  dry  conditions  (Ao).  Other 
parameters  that  were  fixed  included  soil  temperature  (22  oC),  which  defines  the  Henry’s 
Law  Constant,  the  soil  specific  surface  area  (80  m^/g)  which  defines  the  soil  moisture 
content  at  four  monomolecular  layers  (0.089  g/g)  and  the  total  soil  concentration  (Cx= 
1000  ug/kg).  The  chemical  chosen  for  this  simulation  was  2,4-DNT,  as  it  may  likely  be 
one  of  the  best  signature  compounds  for  chemical  detection. 

Two  forecast  evaluations  were  performed.  The  first  was  for  soil  moisture  contents  that 
were  above  the  four  monomolecular  layers  of  soil  water  (0.13  to  0.30  cm?/cm^)  such  that 
the  vapor  solid  partition  coefficient  was  in  the  region  controlled  by  the  HLC.  The 
volumetric  moisture  content  of  0.13  cm?/cm^  is  equivalent  to  a  gravimetric  soil  mositure 
content  of  0.089  g/g  at  a  soil  bulk  density  of  1 .5  g/cnf,  which  is  the  soil  moisture  content 
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at  four  monomolecular  layers  for  a  soil  with  a  specific  area  of  80  m^/g.  The  output 
variable  is  the  soil  gas  concentration.  The  second  forecast  evaluation  was  performed  in 
the  dry  region,  below  four  monomolecular  layers  of  soil  water  (0.01  to  0.13  cmVcm^), 
where  the  soil  vapor  partioning  is  highly  non-linear. 


Assumption:  Soil  Bulk  Density 


Triangular  distribution  with  parameters; 


Minimum 

1.04 

Likeliest 

1.51 

Maximum 

1.68 

Selected  range  is  from  1.04  to  1.68 
Mean  value  in  simulation  was  1.41 

Assumption:  Soil  Particle  Density 

Normal  distribution  with  parameters: 

Mean 

2.60 

Standard  Dev. 

0.01 

Selected  range  is  from  -Infinity  to  +lnfinity 
Mean  value  in  simulation  was  2.60 


Soil  Partide  Density 


327  359  350  353  253 


Assumption:  Soil  Moisture  Content 

Uniform  distribution  with  parameters: 

Minimum  0.13 

Maximum  0.30 


Mean  value  in  simulation  was  0.22 


Soil  Moisture  Content 


0.13  0.17  033  038  050 


Assumption:  Soil  Water  Partition  Coeff  (Kd) 

Lognormal  distribution  with  parameters: 

Mean  3.0 

Standard  Dev.  1.0 

Selected  range  is  from  0.0  to  +lnfinity 
Mean  value  in  simulation  was  3.0 


Soil  Water  Partition  Coeff  (Kd) 


11  2.7  45  55  75 


Assumption:  A  zero 

Uniform  distribution  with  parameters: 

Minimum  10.00 

Maximum  11.00 


A  zero 


lOflO  1035  1050  10.75  11  DO 


Mean  value  in  simulation  was  10.50 
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The  following  forecast  of  soil  gas  concentrations  shows  that  there  is  an  apparent  normal 
or  log  normal  distribution  of  values. 


Figure  5.  Soil  Gas  Concentration  -  Wet 

Figure  6  shows  a  sensitivity  analysis  of  the  input  parameters.  The  strongest  parameter  is 
the  soil  water  partition  coefficient,  which  is  clear  from  a  closer  examination  of  equation 
[9].  The  shape  of  the  histogram  in  Figure  5  is  consistent  with  that  for  the  soil  water 
partition  coefficient,  which  is  also  supported  by  the  strength  of  the  association  as  shown 
in  the  sensitivity  analysis  (Figure  6). 


Sensitivity  Chart 

Target  Forecast  Soii  Gas  Concentration  (ppQ 


Soil  NAaerRatitionODeff  (Md) 

-.99 

Soil  MastuB  Oatent 

.05 

1 

Soii  BJkOeasily 

-.01 

Soil  Particle  Density 

.01 

A2Ero 

.00 

-i  -o'.s  0  0.5  i 

^fe0surBd  by  F^Cortelation 


Figure  6.  Sensitivity  Analysis  -  Wet 

The  next  simulation  was  performed  by  altering  only  the  input  values  for  the  soil  moisture 
content.  In  this  case,  the  soil  moisture  range  is  set  to  be  below  four  monomolecular 
layers  of  water. 
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Assumption:  Soil  Moisture  Content 

Uniform  distribution  with  parameters; 

Minimum  0.01 

Maximum  0.13 


Mean  value  in  simulation  was  0.07 


Soil  Moliture  Content 


Figure  7  shows  that  most  of  the  forecast  values  for  the  soil  gas  concentration  are  well 
below  those  for  the  higher  moisture  content  simulation.  In  addition,  the  sensitivity 
analysis  (Figure  8)  shows  that  the  principal  parameter  is  now  the  moisture  content,  as 
K'd(w)  becomes  dominant  (equations  [9]  and  [12]). 


5,000*61315 


Forecast:  Soil  Gas  Concerrtraf  on  (ppl) 
Frequency  Chart 


52  0udier5 

623 


467.2 


155.7 


0 


Figure  7.  Soil  Gas  Concentration  -  Dry 


Sensitivity  Chart 


Target  Forecast:  Soil  Gas  Concentration  (ppt) 


Soil  Moisture  Content 

.95 

Soil  Water  Partition  Coeff  (Kd) 

-.21 
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Figure  8.  Sensitivity  Analysis  -  Dry 
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This  analysis  shows  that  when  employing  vapor  sniffing  chemical  sensing  technology, 
optimum  conditions  would  be  found  when  the  soils  contained  more  moisture  than  four 
monomolecular  layers. 

B.  1-Dimensional  Transport  Sensitivity  Analysis 

A  one-dimensional  model  sensitivity  analysis  was  used  to  examine  the  impact  of  changes 
to  a  single  parameter  for  a  given  environmental  scenario.  The  output  that  was  evaluated 
was  the  surface  vapor  flux.  The  model  that  was  used  was  the  Buried  Chemical  Model 
from  Jury  et  al.,  1990.  This  sensitivity  analysis  was  presented  in  three  conference 
proceedings  (Appendix  A,  B  and  C)  and  in  a  report  from  Dr.  Jury  (Appendix  D).  Table 
one  shows  the  parameters  evaluated  and  a  summary  of  the  impact. 


Parameter 

Impact  on  steady  state  surface  flux 

Soil  Bulk  Density 

Direct  inversely  proportional 

Henry’s  Law  Constant 

Directly  proportional 

Soil-Water  Partition  Coefficient 

Direct  inversely  proportional 

Source  Flux 

Insignificant  compared  to  initial 
surface  concentration 

Initial  Surface  Concentration 

Directly  proportional 

Burial  Depth 

Increases  lag  time  (very  sensitive) 

Water  Flux  (Precipitation  or  Evaporation) 

Evaporation  enhances,  precipitation 
depresses 

Biochemical  Half-life 

Insensitive  if  >  1  year,  veiy 
sensitive  if  <  60  days 

Dr.  Bill  Jury  performed  a  2-dimensional  analysis,  which  is  included  in  Appendix  E.  This 
showed  that  the  surface  vapor  flux  was  greatest  directly  above  a  source  with  a  small  halo 
up  to  twice  the  length  of  the  buried  source.  However,  the  flux  drops  off  exponentially 
with  lateral  distance. 

C.  Numerical  Model  Development 

A  more  detailed  mechanistic  numerical  model  is  being  developed.  This  model  is  being 
based  on  TOUGH2  (Pruess,  1991)  with  modifications  pertinent  to  the  UXO/landmine 
application  and  is  called  T2TNT.  The  first  round  of  modifications  to  TOUGH2  has  been 
completed  including: 

1.  Addition  of  TNT,  DNT,  and  DNB  vapor  components  -  UXO/landmines  typically 
emit  TNT,  DNT,  and  DNB  vapors.  The  behavior  of  each  of  these  chemicals  is 
different  (vapor  pressure,  vapor/liquid,  liquid/solid,  and  vapor/solid  partitioning),  so 
each  component  is  modeled  separately.  Additional  components  could  be  added  if 
necessary. 

2.  Dusty  Gas  Model  for  gas  diffusion  -  Gas  diffusion  can  be  a  dominant  transport  mode 
for  explosive  vapors  in  the  subsurface,  especially  for  low  moisture  content 
conditions.  In  order  to  mechanistically  model  gas  diffusion  in  a  porous  medium,  the 
Dusty  Gas  Model  (Webb,  1998)  has  been  Implemented. 

3.  Liquid  diffusion  of  dissolved  explosive  gases  -  Liquid  diffusion  can  be  a  dominant 
transport  mode  for  explosive  vapors  in  the  subsurface,  especially  for  moderate  and 
high  moisture  content  conditions.  Liquid  diffusion  was  not  present  in  the  original 
version  of  TOUGH2.  Liquid  diffusion  using  Pick’s  Law  has  been  included  because 
of  the  significant  chemical  concentration  in  the  liquid  phase. 
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4.  Partition  coefficient  as  a  function  of  saturation  -  The  solid  partition  coefficient  may 
be  a  strong  function  of  saturation,  especially  at  low  moisture  content  where  the 
partition  coefficient  may  increase  dramatically  (Petersen,  et  al.  1995).  The  capability 
of  including  a  saturation-dependent  partition  coefficient  has  been  included. 

5.  Biodegradation  -  A  simple  half-life  approach  has  been  implemented  to  model 
biodegradation  of  the  explosive  vapors. 

6.  Surface  Boundary  Conditions  -  Due  to  the  shallow  burial  depth  of  many 
UXO/landmines,  the  fluid  conditions  surrounding  the  UXO/landmine  are  strongly 
influenced  by  the  surface  conditions.  The  parameters  necessary  to  adequately  model 
the  surface  boundary  conditions  are  numerous,  including  solar  and  long-wave 
radiation,  the  surface  boundary  layer  which  is  a  function  of  wind  speed  and  other 
parameters,  precipitation  and  evaporation  at  the  surface,  plants  and  their  root 
systems,  and  the  diurnal  and  seasonal  variation  of  these  parameters.  The  models  used 
for  these  boundary  conditions  are  discussed  in  more  detail  below. 

The  surface  boundary  conditions  discussed  above  are  complex  in  their  own  right. 
Numerous  models  have  been  developed  to  analyze  the  soil-air-plant  system.  In  order  to 
expedite  the  inclusion  of  the  important  surface  conditions  into  T2TNT,  a  number  of 
existing  models  have  been  evaluated.  As  a  result,  the  SiSPAT  model  developed  by  Braud 
et  al.  (1995,  1996)  has  been  selected  for  inclusion  into  T2TNT  with  the  kind  permission 
of  M.  Vauchlin  of  LTHE  in  Grenoble,  France.  Subroutines  from  SiSPAT  have  been 
included  directly  into  T2TNT  as  necessary. 

SiSPAT  has  been  applied  to  a  number  of  field  studies  as  documented  by  Braud  et  al. 
(1995,  1995),  and  Boulet  et  al.  (1997),  and  more  are  in  progress.  Therefore,  SiSPAT 
should  provide  a  well-documented  and  tested  approach  for  modeling  the  soil-plant- 
atmosphere  interface  in  the  T2TNT  code. 

At  the  present  time,  the  surface  boundary  conditions  for  a  bare  soil  have  been 
implemented,  including  the  surface  boundary  layer,  solar  and  long-wave  radiation, 
precipitation,  and  other  conditions  including  the  diurnal  and  seasonal  variation  of  the 
parameters.  Incorporation  of  the  plant  portion  of  the  SiSPAT  model  into  T2TNT  is 
expected  in  1999. 

As  part  of  the  verification  process  for  T2TNT,  comparison  to  results  of  the  one¬ 
dimensional  screening  model  presented  by  Phelan  and  Webb  (1997)  have  been 
performed.  The  conditions  are  for  a  low  desert  environment  with  moderate  moisture 
content;  the  parameters  are  summarized  in  Table  1  as  presented  by  Phelan  and  Webb 
(1997).  The  precipitation/evaporation  cycles  have  not  been  included  in  this  simulation  in 
order  to  provide  comparison  with  the  buried  chemical  analytical  solution  of  Jury  et  al. 
(1990). 

Note  that  some  of  the  assumptions  made  in  the  analytical  model  can  only  be 
approximated  in  the  numerical  code,  such  as  uniform  moisture  content.  In  addition,  there 
are  differences  in  the  gas  diffusion  model,  which  could  lead  to  slightly  different  answers. 
However,  in  general,  the  analytical  and  numerical  problems  are  essentially  equivalent. 
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Table  1.  Simulation  Parameters 


Parameter 

units 

base  case 

e 

cm'*/  cm'* 

0.25 

<|) 

cm^/  cm^ 

0.434 

pb 

g/cm^ 

1.5 

Kd 

cm^/g 

2.5 

Kh 

— 

5.9E-7 

air  boundary 

cm 

0.5 

layer 

hn 

days 

180 

C„ 

pg/cm^ 

4.6E-3 

Jc 

pg/cm^-day 

8.6E-6 

Dr 

cm^/day 

0.432 

Dl 

cmVday 

4320 

burial  depth, 

cm 

5 

top 

burial  depth, 

cm 

15 

bottom 

Figure  9  compares  the  TNT  flux  at  the  surface  from  the  anal3Tical  solution  and  from 
T2TNT.  The  predicted  surface  flux  values  are  higher  than  the  analytical  solution, 
especially  at  smaller  times.  The  difference  is  due  to  numerical  diffusion  of  the  numerical 
method  and  is  expected.  At  longer  times,  the  agreement  between  the  analytical  and 
T2TNT  solutions  is  good. 
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Figure  9.  Comparison  of  TNT  surface  flux  from  Analytical  and  Numerical  Models 
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Figure  10  shows  the  subsurface  distribution  of  TNT  after  1  year  for  the  Jury  et  al.  (1990)  model 
and  from  T2TNT.  The  agreement  in  the  concentration  distribution  is  excellent  at  all  locations. 
Overall,  the  solution  from  T2TNT  agrees  well  with  the  analytical  solution  of  Jury  et  al.  (1990) 
except  for  some  surface  flux  differences  at  early  times  due  to  numerical  diffusion. 

Figure  10.  Comparison  of  Depth  Profile  from  Analytical  and  Numerical  Models 

Additional  verification  and  sensitivity  studies  will  be  performed  with  the  newly- 
developed  model,  including  the  effect  of  the  boundary  layer  on  the  results,  drying 


simulations  including  non-uniform  moisture  content,  diurnal  and  seasonal  variations,  and 
multidimensional  effects  to  determine  the  degree  of  lateral  spreading  including  the 
“footprint”  of  a  landmine. 

Task  Two  —  Fundamental  Property  Measurement 

Early  on  in  this  project,  2,4-DNT  was  identified  as  a  potential  key  chemical  constituent 
that  had  more  favorable  properties  for  chemical  detection,  than  that  of  2,4,6-TNT. 
However,  the  only  water  solubility  vs.  temperature  data  available  was  only  room 
temperature.  We  began  to  evaluate  water  solubility  with  a  method  called  Dynamic 
Coupled  Column  Liquid  Chromatographic  method  that  was  developed  at  the  National 
Bureau  of  Standards  (now  NIST)  for  polycyclic  aromatic  hydrocarbons  in  water  (May 
and  Wasik,  1978).  The  method  is  based  on  generating  saturated  solutions  by  pumping 
water  through  a  column  packed  with  glass  beads  that  have  been  coated  with  the 
compound  of  interest.  This  method  is  recommended  for  low  solubility  compounds  (such 
as  explosives)  and  can  avert  problems  such  as  incomplete  equilibration  of  the  solid  phase 
in  water,  avoidance  of  microcrystals  in  solution  and  adsorptive  loss  on  containers,  filters 
and  transfers  devices.  The  results  from  this  effort  were  not  satisfactory.  The  effluent 
concentration  was  dependent  on  the  flow  rate  passing  through  the  column.  The  results 
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were  not  uniform,  with  some  runs  showing  the  concentration  to  be  a  positive  relationship 
with  flow  rate,  and  others  with  a  negative  slope.  Our  method  used  a  modified  dynamic 
system,  which  used  direct  collection  of  effluent  rather  than  an  extractor  column  as 
described  in  the  published  method.  Subsequent  to  this  effort,  additional  data  for  2,4-DNT 
were  located  (Army,  1971)  and  has  been  used  in  the  modeling  efforts.  However,  this  data 
set  is  lean  and  additional  effort  will  be  spent  to  fill  in  the  data  gaps. 

The  literature  has  shown  the  importance  of  vapor-solid  sorption  for  environmental 
pollutants  such  as  toluene  and  trichloroethene  (Ong  et  al.,  1990).  The  principal  parameter 
for  evaluating  this  phenomenon  is  the  vapor-solid  sorption  coefficient  (K’d).  A  method 
has  been  developed  (Petersen  et  al.,  1988)  that  is  based  on  the  equilibrium  partitioning  in 
closed  spaces  method.  Upon  closer  evaluation  of  this  method  for  compounds  such  as 
2,4,6-TNT  and  2,4-DNT,  it  was  found  that  it  would  be  physically  impossible  due  to  the 
low  vapor  densities.  In  July  1998,  the  DARPA  co-sponsor  of  this  project  requested 
support  for  field  data  collection  of  soil  gas  concentrations  surrounding  buried  landmines. 
The  vapor-solid  sorption  coefficient  data  collection  effort  was  tabled  and  a  change  was 
placed  in  the  September  quarterly  progress  report  pushing  the  milestone  to  December  30, 
1998.  After  consultation  with  the  USA  CRREL,  an  alternate  method  has  been  identified 
that  might  be  suitable  for  this  measurement.  Work  is  currently  in  progress  to  complete 
this  data  collection. 

As  part  of  this  SERDP/DARPA  co-funded  effort,  the  U.S.  Army  Corps  of  Engineers, 

Cold  Regions  Research  and  Engineering  Laboratoiy  has  been  performing  research  on  the 
source  term  from  landmine  materials  and  actual  landmines.  This  effort  has  found  that  the 
initial  surface  concentration  of  explosives  varies  widely,  but  is  about  the  same  magnitude 
as  data  reported  by  Hogan,  1992  of  about  10  ng/cm^.  The  source  flux  values  were 
evaluated  by  two  methods.  One  is  immersion  of  pieces  of  the  plastic  mine  case  in  water 
and  measurement  of  the  concentration  in  water  over  time.  The  other  method  uses  whole 
landmines  in  tedlar  bags  with  measurement  of  the  gas  and  tedlar  bag  sorbed 
concentrations  after  sufficient  length  of  time.  The  values  determined  from  this  effort 
have  shown  source  flux  values  up  to  10^  times  greater  than  that  measured  by  Spangler  et 
al.  1975.  Detailed  technical  reports  on  this  effort  will  be  published  by  Leggett  et  al.  in 
the  future. 

Task  Three  -  Laboratory  Scale  Experiments 

The  laboratory  scale  experiments  have  been  subcontracted  to  New  Mexico  Institute  of 
Mining  and  Technology  (NM  Tech).  This  effort  is  currently  on  schedule  with 
experimental  designs  complete.  The  column  testing  equipment  is  being  acquired  and 
configured  for  use  starting  in  February  1999.  The  basis  for  the  testing  will  be  soil 
columns  with  detailed  soil  moisture  monitoring  with  capability  for  measurement  of 
surface  vapor  flux  and  sorbed  soil  concentrations. 

Task  Four  -  Operational  Strategy 

No  work  was  scheduled  for  this  task  in  1998. 

Task  Five  -  Ordnance  Source  Term 

The  ordnance  source  term  effort  began  in  September  1998.  This  effort  includes 
acquisition  of  live/unfused  ordnance  (105  mm  cartridges,  60  and  81  mm  mortars),  testing 
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for  external  contamination,  performing  immersion  tests  to  evaluate  long  term  source  flux, 
firing  the  ordnance  into  the  ground,  and  performing  post-firing  long-term  source  flux. 

One  pallet  of  each  ordnance  group  (12  individual  ordnance  items)  has  been  located  from 
a  demil  account  at  Redstone  Army  Depot  and  is  currently  in  transit  to  Sandia  National 
Laboratories.  Fuses  will  be  obtained  that  can  be  disassembled  to  remove  the  primary 
charge.  In  this  way,  the  fuse  impact  into  soils  will  be  similar  to  live  UXO. 

The  external  swipe  tests  will  be  performed  in  early  January  1999  and  a  six- week  long 
pre-firing  immersion  test  will  follow.  Careful  consideration  was  taken  in  the  design  of 
the  firing  of  the  cartridges  and  mortars.  Damaged  explosive  devices  have  severe 
constraints  with  transportation  on  public  byways.  Therefore,  collection  of  actual  UXO 
from  current  range  cleanup  operations  and  transportation  to  Sandia  National  Laboratories 
was  not  likely.  However,  Sandia  National  Laboratories  has  the  capabilities  to  fire 
ordnance  into  soils,  recover  the  items,  and  provide  on-site  transportation. 

Two  designs  were  evaluated  for  firing  the  ordnance  into  soils.  An  above  ground  soil 
berm  of  compacted  soil  with  concrete  containment  walls  was  considered;  however, 
consultation  with  penetrator  testing  groups  at  Sandia  National  Laboratories  discouraged 
this  method  as  containment  is  typically  not  guaranteed.  Using  an  Army  manual 
(Technical  Manual  TM  5-855-1,  July  1965),  penetration  of  the  105  mm  cartridge  is 
expected  to  be  about  7  to  9  feet  at  an  impact  velocity  of  800  to  900  fps.  A  low  angle  (15 
to  20°)  shot  into  native,  undisturbed  soil  with  a  penetration  depth  of  20  feet  yields  a  depth 
below  grade  of  about  10  feet.  By  probing  the  soil  trajectory  path,  the  angle  and  total 
penetration  length  can  be  determined,  and  recovery  of  the  target  can  be  performed  with 
traditional  excavation  equipment.  Explosive  Ordnance  Disposal  trained  services  area 
available  to  recover  the  fired  ordnance.  Therefore,  a  trench  will  be  excavated  that  will 
serve  as  the  face  for  the  low  angle  shot  into  undisturbed  soil. 
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List  of  Symbols 

Pb-  soil  bulk  density  (g/cm^) 

Ps  -  soil  particle  density  (g/cm^) 

Pw  -  density  of  water  (1.0  g/cm^) 

0  -  volumetric  moisture  content  (cmVcm^) 
w  -  gravimetric  soil  moisture  content  (g/g) 

W4  -  water  content  at  4  monomolecular  layers  of  water  (g/g) 

<p  -  soil  porosity  (cmVcm^) 
a  -  air  filled  porosity  (cmVcm^) 

Y  -  activity  of  water  (~1) 

Ct  -  Total  soil  concentration  (g/cm^) 

Cs  -  Concentration  on  solid  phase  (g/cm^) 

Cl  -  Concentration  in  soil  aqueous  phase  (g/cm^) 

Cg-  Concentration  in  the  soil  gas  (g/cm^) 

CsG  -  Concentration  on  soil  solid  phase  from  vapor  sorption  (g/cm^) 
CsL  -  Concentration  on  soil  solid  phase  from  liquid  sorption  (g/cm^) 
Rs  -  Soil  solid  phase  partitioning  factor  (g/cm^) 

Rg  -  Soil  gas  phase  partitioning  factor  (cm^/cm^) 

Rl  —  Soil  liquid  phase  partitioning  factor  (cmVcm^) 

Rsg-  Soil  solid-vapor  phase  partitioning  factor  (g/cm^) 

Rsl  -  Soil  solid-liquid  phase  partitioning  factor  (g/cm^) 


Kd—  soil  water  partition  coefficient  (cmVg) 

Kh  -  Henry’s  Law  Constant  (cm^air/cm^water) 

Ksg  —  soil  vapor  partition  coefficient  (cmVg) 

K'd  -  soil  vapor  partition  coefficient  (cmV|) 

Ao  -  log  of  K'd  at  oven  dry  conditions  (cmvg) 

A4  -  log  of  K'd  at  moisture  content  of  4  monomolecular  layers  of  water 
P(w)  —  fitting  parameter 
a  -  fitting  parameter  for  K'd(w) 

Sa  -  soil  specific  surface  area  (m^/g) 

MWw  -  molecular  weight  of  water  (18  g/mole) 

MA„  —  molecular  area  of  water  (10.8E-20  mVmolecule) 

An  —  Avogadro’s  Number  (6.02E23  molec/mole) 
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Chemical  Detection  of  Buried  Landmines 


James  M.  Phelan*  and  Stephen  W.  Webb^ 


Of  all  the  buried  landmine  identification  technologies  currently 
available,  sensing  the  chemical  signature  from  the  explosive 
components  found  in  landmines  is  the  only  technique  that  can 
classify  non-explosive  objects  from  the  real  threat  In  the  last  two 
decades,  advances  in  chemical  detection  methods  have  brought 
chemical  sensing  technology  to  the  foreground  as  an  emerging 
technological  solution.  In  addition,  advances  have  been  made  in 
the  understanding  of  the  fundamental  transport  processes  that 
allow  the  chemical  signature  to  migrate  from  the  buried  source  to 
the  ground  surface.  A  systematic  evaluation  of  the  transport  of 
the  chemical  signature  from  inside  the  mine  into  the  soil 
environment,  and  through  the  soil  to  the  ground  surface  is  being 
explored  to  determine  the  constraints  on  the  use  of  chemical 
sensing  technology.  This  effort  reports  on  the  results  of 
simulation  modeling  using  a  one-dimensional  screening  model  to 
evaluate  the  impacts  on  the  transport  of  the  chemical  signature  by 
variation  of  some  of  the  principal  soil  transport  parameters. 

landmines,  chemical  sensors,  soil  transport 

L  Introduction 


The  organic  chemicals  of  the  explosives  in  the  buried 
landmine  environment  can  exist  in  or  on  four  phases:  solid 
phase  of  the  neat  explosive  material,  vapor  phase  in  the  soil 
air,  aqueous  phase  in  the  soil  w'ater  solution,  and  sorbed 
onto  soil  solid  phases.  The  chemical  signature  begins  as  a 
surface  coating  from  production  or  depot  storage  and 
through  continuous  emission  by  permeation  through  the 
mine  case  or  through  leaks  in  seals  and  seams.  Once  the 
chemicals  enter  the  soil  environment,  they  experience  phase 
transitions,  partitioning  into  the  soil  air,  soil  water  and 
sorbing  onto  soil  particles.  The  impact  of  temperature  and 
chemical  gradients,  and  precipitation/evaporation  will  cause 
movement  of  the  chemical  signature.  Part  of  this  transport 
is  upward  to  the  soil  surface  where  chemical  detection 
technology  is  envisioned  to  be  used.  Simulation  modeling 
is  a  technique  that  can  evaluate  the  impacts  of  many  of  the 
environmental  variables  that  can  dampen  or  accentuate  the 
surface  expression  of  the  chemical  signature.  Model  results 
will  be  shown  that  describe  the  magnitude  of  the  changes 
that  accompany  variations  due  to  chemical  properties  of  the 
explosive  and  properties  of  the  soil  environment. 


*  Environmental  Restoration  Technologies  Department, 
^Mission  Analysis  and  Simulation  Department, 
at  Sandia  National  Laboratories,  Albuquerque,  NM. 


11.  Chemicals  in  the  Soil  environment 

Soils  are  porous  media  with  a  number  of  physico-chemical 
properties  that  affect  the  transport  of  explosive  chemicals. 
Soil  bulk  density  is  a  measure  of  the  compaction  of  the  soil 
and  is  defined  as 

M, 


Pi=- 


[1] 


where  pj,  is  the  soil  bulk  density  (g/cm^),  M,  is  the  mass  of 
soil  particles  (g),  and  V,  is  the  volume  of  soil  (cm^).  Soils 
under  natural  conditions  have  bulk  densities  ranging  from 
1.0  to  1.8  g/cm^.  However,  soils  that  have  been  excavated 
and  replaced,  such  as  during  the  emplacement  of  a 
landmine,  may  have  bulk  densities  much  less  than  1,  The 
soil  bulk  density  is  inversely  proportional  to  the  soil 
porosity  as  follows 

(j>  =  l-p,/p,  [2] 

where  p,is  the  soil  particle  density  (ranges  from  2.6  to  2.8 
g/cm’  for  most  soils).  The  soil  porosity,  or  void  volume,  is 
defined  as 

V  +v 

[3] 

"  s 

where  <j)  is  the  soil  porosity  (emVem^),  is  the  volume  of 
soil  water  (cm^)  and  V,  is  the  volume  of  soil  air  (cm^).  Soil 
porosity  values  range  from  0.3  for  sands  to  0.6  for  clay  rich 
soils.  The  volumetric  moisture  content  describes  how 
much  water  is  present  in  the  soil  and  changes  greatly  during 
precipitation/drainage  events  and  evaporation  conditions. 
Volumetric  water  content  is  defined  as 

S 

where  6  is  the  volumetric  water  content  (emVem^).  Soil 
moisture  contents  have  values  from  near  zero  up  to  the  soil 
porosity  value.  When  the  soils  are  not  fully  saturated,  the 
balance  of  the  soil  pore  space  not  filled  with  water  is 
termed  the  air  filled  porosity,  and  is  defined  as 


a  = 


[5] 


where  is  the  volumetric  air  content  (emVem^). 


It  is  often  more  convenient  to  use  soil  saturation  (S^) 
because  it  is  a  measure  of  the  relative  saturation  of  a 
particular  soil  pore  space  with  water. 


e 

Since  the  explosive  chemicals  can  exist  as  solutes  in  the  soil 
water  and  the  movement  of  soil  water  can  be  a  significant 
transport  mechanism,  water  solubility  is  an  important 
parameter.  Water  solubility  is  defined  as 


where  Q  is  the  concentration  in  aqueous  phase 
(g/cm^  soil  water)  and  is  the  mass  of  chemical  (e.g. 
TNT)  (g).  Water  solubility,  however,  is  not  constant  and  is 
typically  an  increasing  flmction  with  temperature. 


Henry’s  Law  constant  is  a  relative  measure  of  the  amount 
of  the  chemical  that  exists  in  the  gas  phase  to  that  in  the 
aqueous  phase  at  equilibrium,  and  is  defined  as 

[8] 

where  Kh  is  the  Henry’s  Law  constant  (unitless)  and  Co  is 
the  concentration  in  gas  phase  (g/cm3  soil  gas).  Henry’s 
Law  constant  is  also  a  function  of  temperature  because  both 
Cq  and  Cl  are  functions  of  temperature. 


The  soil  partition  coefficient  is  a  relative  measure  of  how 
much  of  the  chemical  is  temporarily  bound  to  the  soil  to 
that  in  the  soil  aqueous  phase 

^.=7f  [9] 

where  IQ  is  the  linear  soil-water  partition  coefficient 
(cmVg)  for  water  saturated  soils  and  Cg  is  the  concentration 
sorbed  on  the  soil  solid  phase  (g/g  of  soil). 


The  soil  water  partition  coefficient  is  often  correlated  with 
the  fraction  of  organic  carbon  found  in  the  soils.  In  this 
way,  the  variability  between  soils  can  be  reduced.  The 
organic  carbon  distribution  coefficient  is  defined  as 

K^=-f  [10] 

J  oc 

where  is  the  organic  carbon  distribution  coefficient  and 
is  the  fraction  of  organic  carbon. 


Table  1.  Vapor  Density,  Water  Solubility  and  Henry’s  Law 


Constant  of  Explosive  Compounds  at  20'^C 


TNT 

DNT 

RDX 

Vapor  Density  (pg/m^) 

43.5 

122 

0.024 

Water  Solubility  (mg/1) 

130 

270 

50 

Henry’s  Law  Constant 

3.35E-7 

4.51E-7 

4.73E-10 

Pennington  and  Patrick  (1990)  measured  the  soil  water 
partition  coefficient  (IQ)  of  TNT  in  fourteen  soils  from 
locations  across  North  America.  The  mean  value  was  3.8 
cmVg  with  a  standard  deviation  of  1.34.  The  highest  value 
was  6.8  cmVg  and  the  lowest  value  was  2.3  cmVg.  Xue  et 
al.  (1995)  evaluated  two  soils  and  showed  mean  values  for 
TNT  of  2.66  cmVg  and  3,64  cmVg.  DNT  and  RDX  have 
very  little  data.  Xue  et  al.  (1995)  showed  values  for  RDX 
of  1.59  cmVg  and  1.57  cmVg.  McGrath  (1995)  showed  a 
IQ,  value  of  25 1  for  DNT.  For  the  fourteen  soils  evaluated 
by  Pennington  and  Patrick  (1990),  the  mean  value  for  the 
fraction  of  organic  carbon  was  0.0173  with  a  standard 
deviation  of  0.01 1.  Using  these  values,  the  IQ  for  DNT  has 
a  mean  value  of  4.4±2.7  cmVg  (one  std.  dev.).  In  summary, 
the  soil  water  partitioning  coefficients  for  TNT,  DNT  and 
RDX  all  fall  into  an  approximate  range  between  1,5  and  7.0 
cmVg.  This  is  a  rather  narrow  range  as  common  chemicals 
can  have  values  one  to  two  orders  of  magnitude  lesser  and 
greater  than  these. 

The  biochemical  half-life  of  explosives  in  near  surface  soils 
has  not  been  studied  well  outside  of  the  biotreatment 
technology  area  for  contaminated  soils.  However,  long¬ 
term  surface  soil  degradation  tests  at  Los  Alamos  National 
Laboratory  (LANL)  followed  the  degradation  of  soils 
doped  with  1000  mg/kg  of  various  explosives  over  20  years 
(Dubois  and  Baytos,  1991).  Table  2  shows  the  half-lives 
estimated  from  these  long-term  experiments. 


Table  2.  Estimates  of  Half-Lives  of  Explosives  from 
LANL  Long-Term  Surface  Soil  Tests _ 


Explosive 

Half-Life  (years) 

TNT 

1 

RDX 

36 

HMX 

39 

PETN 

92 

III.  Physical  Chemical  Properties  of  Explosives 

The  principal  explosive  chemicals  found  in  landmines  are 
TNT  and  RDX  (NGIC,  1995).  DNT,  as  a  production  by¬ 
product  of  TNT,  is  also  considered  to  be  a  significant 
signature  chemical  for  buried  landmines.  As  a  group,  these 
chemicals  have  very  low  vapor  densities  and  moderately 
low  water  solubilities.  Table  1  shows  these  properties  and 
the  Henry’s  Law  constant  at  20®C  (Phelan  and  Webb, 
1997). 


IV.  Screening  Model 

The  environmental  fate  and  transport  of  organic  chemicals 
including  volatilization  and  leaching  losses  has  been  used 
to  explore  the  distribution  of  agricultural  pesticides  in  soils 
(Mayer  et  al.  1974,  Farmer  et  al,  1980,  and  Jury  et  al. 
1980).  These  models  were  primarily  intended  to  simulate 
specific  circumstances.  However,  Jury  et  al.  (1983, 1984a, 
1984b,  1984c)  developed  and  validated  a  general  screening 
model  (Behavior  Assessment  Model,  BAM)  that  included 
volatilization,  leaching,  and  degradation  to  explore  the 
major  loss  pathways  of  agricultural  pesticides  as  a  function 


of  specific  environmental  conditions.  The  model 
simulations  can  be  used  to  assess  the  behavior  of  different 
chemicals  under  particular  environmental  conditions,  but  is 
not  intended  to  predict  a  definitive  concentration 
distribution  in  the  field.  As  such,  the  predictions  from  the 
screening  model  are  only  an  indication  of  expected 
conditions. 

This  model  is  valuable  in  that  it  can  express  the  total 
concentration  of  a  chemical  in  the  gas,  aqueous  and  sorbed 
phases.  The  total  concentration  is  expressed  as 


4.  The  chemical  will  undergo  first-order  degradation  due 

to  biological  and  chemical  effects. 

5.  Chemical  movement  is  one  dimensional. 

6.  The  adsorbed  and  dissolved  phases  undergo 

reversible,  linear  adsorption. 

7.  The  dissolved  and  gaseous  phase  concentrations  are 

related  through  Henry’s  law. 

8.  The  soil  properties  are  constant  in  space  and  time. 

9.  Water  flux  is  constant  in  space  and  time  (relaxed  in 

the  present  application). 

10.  Volatilization  of  the  chemical  to  the  atmosphere  is  by 
vapor  diffusion  through  an  air  boundary  layer  of 
constant  thickness. 


where  Cg  is  the  concentration  sorbed  to  the  soil.  Cl  is  the 
solute  concentration  in  the  aqueous  phase,  and  Cq  is  the  gas 
phase  concentration.  In  addition.  Jury  (1983)  shows  how 
equation  [1 1]  can  be  rewritten  in  terms  of  one  of  the 
variables  alone 


[12] 


where 


„  e 

=  Pa  +  ^ 


Em 

K, 


[13] 


Rj^  =  Pi,K^  +9  +  5 


=  Pa 


K,  e 


[14] 

[15] 


are  the  solid,  liquid  and  gas  phase  partition  coefficients, 
respectively. 


An  adaptation  of  the  BAM  was  developed  to  be  applicable 
to  the  conditions  of  contaminated  soil  buried  under  a 
known  depth  of  clean  soil  -  Buried  Chemical  Model,  BCM 
(Jury  et  al.,  1990).  Simulations  based  on  a  modification  of 
Jury’s  BCM  are  used  in  this  report  to  simulate  the  behavior 
of  the  chemical  signature  from  buried  landmines.  The 
Buried  Chemical  Model  of  Jury  et  al.  (1990)  is  based  on  the 
following  assumptions.  A  detailed  discussion  of  these 
assumptions  is  given  in  Jury  et  al  (1990). 


1.  The  chemical  may  adsorb  on  the  solid  phase,  be 

dissolved  in  the  aqueous  phase,  or  exist  in  the  vapor 
phase. 

2.  The  chemical  flux  is  the  sum  of  the  vapor  flux  and  the 

dissolved  solute  flux  using  Pick’s  law. 

3.  The  porous  medium  factors  for  gas  and  liquid  phase 

diffusion  are  given  by  the  Millington  and  Quirk 
(1961)  model  as  extended  for  liquid  diffusion  by 
Jury  et  al.  (1983). 


In  the  present  implementation  of  Jury’s  model,  a  constant 
source  term  has  been  added  to  reflect  the  chemical  source 
from  the  landmine  at  a  specific  location. 


Under  these  assumptions  (including  the  source  term)  the 
model  formulation  becomes 


dCr 


+  \iCf  — 


d^Cr 


dCr 


[16] 


dt  ■  dz^  "  dz 

where  Cr  is  the  total  chemical  concentration,  \x  is  the  bio¬ 
chemical  decay  constant,  and  o  is  the  source  term.  The 
effective  velocity  (V^)  is  defined  as 

f [17] 

Pb^d  +0 

where  J^,  is  the  precipitation/evaporation  flux.  The 
effective  diffusion  coefficient  (Dg)  of  the  chemical  is 
defined  as 


,10/3  i 


klO/3 


d; 


[18] 


where  D°  is  the  diffusivity  of  the  gas  phase  of  the 


chemical  in  air  and  is  the  diffusivity  of  the  chemical  in 
aqueous  phase .  The  boundary  conditions  for  the  problem 
are  diffusion  through  a  boundary  layer  at  the  upper  surface, 
and  a  zero  chemical  concentration  at  infinity  at  the  lower 
boundary.  These  boundary  conditions  can  be  expressed  as 


-A 


ac 


dz 


^-+V,Cr  =  -H,C, 


where 


and 


hK 


H 


h  = 


d; 


[19] 


[20] 


[21] 


and 


=  0- 


The  initial  conditions  are  an  initial  concentration,  Cq,  over 
an  interval  from  L  to  W,  or 

Q(z,0)  =  0  L)z)fV 

Q(2,0)  =  Co  L<z<W 

The  above  model  results  in  a  closed  form  solution  as  a 
function  of  space  and  time;  the  results  are  rather  lengthy 
and  will  not  be  presented  here  but  are  given  by  Jury  et  al. 
(1983, 1990).  In  the  present  simulation,  the  assumption  of 
constant  water  flux  in  time  will  be  relaxed.  Therefore, 
sequences  of  water  fluxes  representing  desired  conditions 
(rainfall  followed  by  evaporation)  can  be  simulated  to 
determine  the  effect  of  water  flux  variations  on  the  location 
of  TNT  in  the  soil  and  the  surface  TNT  vapor  flux.  A 
numerical  solution  was  developed  and  verified  by 
comparison  to  the  results  given  by  Spencer  et  al.  (1988) 
and  Jury  et  al.  (1990)  (Phelan  and  Webb,  1997). 

Using  this  solution,  simulations  were  performed  using  a 
landmine  that  has  contributed  an  initial  soil  concentration 
(Cp)  based  on  the  surface  contamination  of  the  landmine.  It 
has  been  assumed  that  the  entire  surface  contamination  was 
completely  and  uniformly  transferred  to  the  soil  just  prior 
to  the  beginning  of  the  simulation  runs.  Surface 
contamination  data  (Hogan  et  al.,  1992)  showed  a  median 
surface  contamination  of  15  ng/cm^  from  42  domestic  and 
foreign  landmines.  Using  the  dimensions  of  an  anti-tank 
(AT)  mine  of  30  cm  diameter  by  10  cm  high,  the  surface 
contamination  would  provide  3.5x10'^  g  of  TNT  for  initial 
distribution  in  the  soil.  Using  the  volume  of  the  AT  mine 
that  this  mass  of  TNT  is  distributed  into,  the  initial 
concentration  (CJ  would  be  -5x10'^  |ag/cm^ 

The  constant  source  term  emanation  rate  was  derived  from 
vapor  collection  chamber  experiments  on  two  mines 
(Spangler,  1975).  Values  ranged  from  10‘*^to  10‘*®g/cm^- 
s.  The  higher  rate  of  10’*^  g/cm^-s  (8.6x10'^  pg/cm^-day) 
was  used  in  these  simulations.  If  the  top  of  the  AT  mine 
was  buried  at  a  depth  of  10  cm,  the  burial  zone  of  the  initial 
contamination  is  from  10  to  20  cm,  and  the  constant  source 
term  is  placed  at  a  depth  of  15  cm. 

The  diffusivity  of  gas  in  air  (  )  and  diffusivity  of  liquid 

in  water  (  )  were  selected  from  Jury  et  al.  (1983).  The 

biochemical  half-life  value  of  365  days  was  selected  from  a 
long  term  field  experiment  (Dubois  and  Bayton,  1991). 

The  precipitation/evaporation  rates  and  periods  followed  in 
all  the  simulations  here  were  the  low  desert  scenario  from 
Phelan  and  Webb  (1997).  This  scenario  was  derived  from 
data  found  in  HELP  (Hydrological  Evaluation  of  Landfill 
Performance)  model  (Schroeder  et  al.,  1994a  and  1994b). . 
The  HELP  model  showed  that  the  low  desert  had  1  day  of 
precipitation  followed  by  7  days  of  evaporation.  For 


simplicity,  total  precipitation  and  total  evaporation  for  each 
cycle  are  assumed  to  be  equal  and  for  these  simulations  the 
cycles  were  continued  for  approximately  four  to  ten  years. 

Table  2  shows  the  input  parameters  used  in  the  simulations. 


Table  2.  Simulation  Parameters 


parameter 

units 

base  case 

variant 

cases 

e 

cmV  cm^ 

0.25 

0.375 

<l> 

cmV  cm’ 

0.5 

♦ 

pb 

g/cm^ 

1.5 

♦ 

K, 

cmVg 

1.6 

3.8 

6.0 

Kh 

5.9E-7 

4.73E-10 

air  boundary 

cm 

0.5 

♦ 

layer 

h/2 

days 

365 

180 

60 

c. 

pg/cm^ 

4.6E-3 

0 

Jc 

pg/cm^-day 

8.6E-6 

0 

d: 

cmVday 

0.432 

♦ 

cmVday 

4320 

♦ 

burial  depth, 

cm 

10 

♦ 

top 

burial  depth. 

cm 

20 

♦ 

bottom 

cycles/yr 

- 

45 

* 

precipitation 

days 

1 

* 

evaporation 

days 

7 

* 

precipitation 

cm/day 

0.44 

* 

rate 

evaporation 

cm/day 

-  0.063 

♦ 

rate 

total 

cm/year 

20 

precipitation/ 

evaporation 

*  -  same  as  the  base  case 


V.  DISCUSSION 

For  each  of  the  figures  shown,  there  is  a  distinct  oscillation 
of  the  surface  vapor  flux.  This  feature  is  a  result  of  the 
cycling  of  precipitation  and  evaporation.  To  evaluate  the 
effect  of  the  Henry’s  Law  constant,  two  simulations  were 
performed  where  all  parameters  were  kept  constant  with 
one  case  using  a  Kh  equal  to  that  of  TNT  and  one  for  RDX 
(both  at  20°C).  Figure  1  shows  the  results  and  indicate  that 
the  TNT  surface  vapor  flux  would  be  expected  to  reach  a 
greater  steady  state  value  than  RDX,  approximately 
proportional  to  the  ratio  of  Henry’s  Law  constants.  For 
TNT,  a  temperature  increase  from  0°C  to  40°C  will 
increase  the  value  by  a  factor  of  about  100.  It  appears 
that  seasonal  and  diurnal  soil  temperature  changes  could 


make  a  significant  effect  on  the  subsurface  transport  and 


surface  flux  of  explosive  signatures. 


Tim#  (days) 


Figure  1.  Effect  of  Henry’s  Law  Constant  on  Surface  Vapor  Flux 

Depth  concentration  profiles  at  the  end  of  the  simulation  Next,  simulations  were  performed  to  evaluate  small 

period  show  that  for  both  cases  the  concentration  of  TNT  changes  in  the  soil  water  partitioning  coefficient.  The 

and  RDX  are  essentially  equal.  This  implies  that  most  of  values  for  TNT,  DNT  and  RDX  have  values  from  about  1.5 

the  transport  upward  to  the  ground  surface  is  within  the  to  7.  Figure  2  shows  the  surface  flux  over  time  for 

aqueous  phase  and  that  the  release  of  the  chemical  into  the  values  of  1.6,  3.8  and  6.0  cmVg. 

vapor  phase  above  the  ground  surface  is  directly 
proportional  to  the  Henry’s  Law  constant. 


Tmw  (days) 

Figure  2.  Effect  of  Soil  Water  Partitioning  Coefficient  on  Surface  Vapor  Flux 

These  simulations  show  that  even  though  the  soil  water  steady  state  concentrations  stabilize  at  much  lower  levels, 

partitioning  coefficients  appear  to  vary  only  slightly  among  Also  with  the  lower  Kj  values,  the  effect  of  precipitation/ 

many  soils,  there  is  a  significant  impact  to  the  transport  of  evaporation  cycles  becomes  more  pronounced.  This  is 

the  chemical  to  the  ground  surface.  As  the  K,,  value  consistent  with  the  lower  value,  since  more  of  the  mass 

increases,  the  lag  period  becomes  much  longer  and  the  of  the  chemical  is  found  in  the  aqueous  phase  and  is 


affected  by  the  upward  and  downward  flux  of  water. 
Figure  3  shows  the  depth  profiles  from  the  simulations 
varying  the  soil  water  partitioning  coefficients.  These 


curves  show  that  the  simulations  with  lower  values  have 
more  significant  transport  of  the  chemicals  to  soils  both 
above  and  below  the  source  zone  than  the  higher  values. 


Figure  3.  Depth  Profile  After  4  Years  With  Variant  Values 


The  next  simulations  were  an  evaluation  of  the  source  term 
parameters.  Very  little  data  is  available  on  the  initial 
concentration  of  explosives  coating  the  outer  surface  of  a 
mine  and  even  less  so  on  the  surface  emission  flux.  Figure 
4,  curve  one  shows  the  effect  of  reducing  the  surface 
emission  flux  to  essentially  zero  (a  value  of  lE-20  ug/cm2- 
day  was  used).  This  curve  is  essentially  the  same  as  curve 
one  in  Figure  1  where  the  surface  emission  flux  (1^)  was 
equal  to  8,6E-6  ug/cm2*day.  This  implies  that  the  surface 
emission  flux  makes  very  little  contribution  to  the  overall 
mass  transport  in  the  soil.  If  one  considers  the  total  mass 
of  explosive  contributed  by  the  initial  concentration  (35  ug) 


and  the  amount  from  the  surface  emission  flux  (0.02 
ug/day),  it  would  take  about  4.8  years  for  the  total  mass  of 
explosive  from  the  surface  emission  flux  to  equal  the 
amount  from  the  initial  surface  contamination.  Curve  two 
in  Figure  4  shows  the  effect  where  the  initial  concentration 
on  the  mine  is  essentially  zero  (a  value  of  lE-20  ug/cm3 
was  used).  This  shows  a  significantly  longer  lag  time  and 
about  four  orders  of  magnitude  lower  steady  state  surface 
vapor  flux  at  the  end  of  the  simulation  period.  Figure  4 
indicates  that  the  initial  concentration  is  a  much  more 
important  parameter  than  the  surface  flux  for  the  mass 
transport  of  chemicals  to  the  ground  surface. 


Figure  4.  Effect  of  Source  Term  Variation  on  Surface  Vapor  Flux 


Figure  5  shows  the  cases  where  the  initial  concentration  is 
increased  by  a  factor  of  1 0  (curve  2)  and  a  factor  of  1 OE+5 
(curve  1)  over  the  base  case  (curve  3).  These  simulations 


show  that  the  lag  period  remains  about  the  same;  however, 
the  increase  in  steady  state  flux  is  proportional  to  the 
increase  in  initial  concentration  (Q). 
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Figure  5.  Effect  of  Source  Term  Variation  on  Surface  Vapor  Flux 

Another  parameter  where  there  is  very  little  data  is  the  (curve  2),  Figure  6  shows  that  over  the  long-term,  the 

biochemical  half  life  (tj/^).  There  are  many  influences  on  shorter  half-life  will  significantly  decrease  the  steady  state 

the  magnitude  of  this  parameter  and  the  variability  is  surface  flux.  Biochemical  decay  constants  that  are  very 

expected  to  be  large.  Simulations  over  ten  (10)  years  were  large  (e.g.  RDX  values  over  30  years)  appears  to  have 

completed  to  assess  the  impact  of  decreases  in  the  minimal  impact  to  the  short  term  soil  transport  phenomena 

biochemical  half-life  from  365  days  (curve  1)  to  180  days  evaluated  in  these  simulations. 
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Figure  6.  Effect  of  Half-Life  on  Surface  Flux 

Initial  simulations  using  this  approach  showed  that  for  a  concentrations  (Phelan  and  Webb,  1997).  Higher  moisture 

particular  scenario,  the  constant  moisture  content  value  had  contents  showed  significantly  shorter  lag  periods  and 

a  significant  impact  on  the  lag  period  (when  the  vapor  flux  greater  steady  state  concentrations.  Figure  7  shows  how 

reached  the  ground  surface)  and  the  steady  state  increasing  the  soil  saturation  from  0.5  to  0.75,  decreases  the 


lag  period  substantially  and  increases  the  steady  state  simulations.  Future  efforts  will  include  the  development  of 

surface  flux  by  two  orders  of  magnitude.  The  constant  a  numerical  simulation  capability  that  can  explore  the 

moisture  content  assumption  in  this  screening  model  allows  transport  phenomena  of  chemical  signatures  from  buried 

for  a  simplified  mathematical  solution  to  complex  transport  landmines  in  more  detail  and  with  fewer  simplifying 

phenomena.  One  must  recognize  this  assumption  and  not  assumptions, 

over  generalize  the  information  gleaned  from  these 
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Figure  7.  Effect  of  Soil  Moisture  Content  on  Surface  Flux 

In  some  parts  of  the  world,  minefields  are  located  in  areas  followed  by  -0.5  cm/day  of  evaporation  for  60  days.  This 

that  experience  very  distinct  wet  and  dry  climatic  periods.  shows  the  immediate  drop  in  the  steady-state  surface  flux 

Figure  8  shows  the  result  of  a  simulation  to  assess  the  effect  after  precipitation  begins.  Once  the  evaporation  period 

of  a  short-term  continuos  precipitation  period  followed  by  a  begins,  there  is  a  short  lag  period  where  the  surface  vapor 

short-term  evaporation  period.  The  baseline  simulation  flux  stays  nearly  constant  before  rising  to  just  above  the 

(Figure  1,  curve  1)  was  run  for  1440  days  followed  by  1  flux  before  the  precipitation  began, 

cm/day  of  precipitation  for  30  days  which  was  then 
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Figure  8.  Effect  of  Continuous  Precipitation  Followed  by  Continuous  Evaporation 


VI.  Summary 

The  detection  of  buried  landmines  has  become  a  significant 
challenge  to  the  technical  community.  The  challenges  of 
locating  a  buried  object  in  near  surface  soils  have  driven  the 
need  for  technology  from  traditional  geophysical  sensor 
systems  to  those  looking  for  the  chemical  signatures  of  the 
explosives  derived  from  within  the  landmine.  The  transport 
of  the  chemicals  found  in  the  explosive  charge  through 
soils  is  a  complex  process  involving  phase  changes, 
interactions  with  the  soils,  and  biochemical  reactions. 

An  investigation  with  computational  simulation  was 
initiated  to  explore  the  impacts  of  several  of  the  various 
input  parameters  with  a  pesticide  screening  model  adapted 
to  the  landmine  chemical  sensing  problem.  It  was  found 
that  the  variations  in  the  Henry’s  Law  constant  are  directly 
proportional  to  the  changes  seen  in  the  steady-state  surface 
vapor  flux.  Small  changes  in  the  soil  water  partitioning 
coefficient  made  large  changes  in  the  lag  period  and  steady- 
state  surface  vapor  flux.  The  initial  source  term  and 
continuous  source  flux  of  chemicals  from  the  landmines  is 
poorly  understood  and  expected  to  be  highly  variable. 

These  simulations  found  that  with  the  baseline  case,  if  the 
continuous  source  flux  was  absent,  there  was  no  significant 
difference  in  the  surface  vapor  flux  or  subsurface 
distribution  at  the  end  of  the  simulation  period.  This 
implies  that  the  continuous  source  flux  may  be  much  less 
important  than  the  initial  surface  contamination.  It  appears 
that  the  magnitude  of  the  surface  vapor  flux  is  directly 
proportional  to  the  amount  of  the  initial  surface 
contamination.  The  biochemical  half-life  is  another 
parameter  that  is  likely  to  have  very  different  values 
depending  on  the  location  and  climatic  conditions. 
Simulations  over  ten  years  showed  that  the  steady  state 
surface  vapor  flux  declines  steadily  when  the  biochemical 
half-life  becomes  smaller  than  one  year.  Finally,  the 
importance  of  heavy  precipitation  (such  as  a  monsoon 
season)  followed  by  a  dry  season  was  explored.  The 
impact  of  the  heavy  precipitation  was  to  lower  the  surface 
flux  seven  orders  of  magnitude;  however,  the  evaporation 
period  that  followed  returned  the  surface  flux  to 
approximately  the  pre-monsoon  surface  flux. 
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ABSTRACT 

The  fate  and  transport  of  chemical  signature  molecules  that  emanate  from  buried  landmines  is  strongly  influenced  by 
physical  chemical  properties  and  by  environmental  conditions  of  the  specific  chemical  compounds.  Published  data  have 
been  evaluated  as  the  input  parameters  that  are  used  in  the  simulation  of  the  fate  and  transport  processes.  A  one-dimensional 
model  developed  for  screening  agricultural  pesticides  was  modified  and  used  to  simulate  the  appearance  of  a  surface  flux 
above  a  buried  landmine  and  estimate  the  subsurface  total  concentration.  The  physical  chemical  properties  of  TNT  cause  a 
majority  of  the  mass  released  to  the  soil  system  to  be  bound  to  the  solid  phase  soil  particles.  The  majority  of  the  transport 
occurs  in  the  liquid  phase  with  diffusion  and  evaporation  driven  advection  of  soil  water  as  the  primary  mechanisms  for  the 
flux  to  the  ground  surface.  The  simulations  provided  herein  should  only  be  used  for  initial  conceptual  designs  of  chemical 
pre-concentration  subsystems  or  complete  detection  systems.  The  physical  processes  modeled  required  necessary 
simplifying  assumptions  to  allow  for  analytical  solutions.  Emerging  numerical  simulation  tools  will  soon  be  available  that 
should  provide  more  realistic  estimates  that  can  be  used  to  predict  the  success  of  landmine  chemical  detection  surveys  based 
on  knowledge  of  the  chemical  and  soil  properties,  and  environmental  conditions  where  the  mines  are  buried.  Additional 
measurements  of  the  chemical  properties  in  soils  are  also  needed  before  a  fully  predictive  approach  can  be  confidently 
applied. 

Keywords:  Mine  detection,  chemical  sensing,  soil  transport,  computational  simulation 

1.  INTRODUCTION 

The  goal  of  locating  buried  landmines  is  a  significant  challenge  to  science  and  technology  (Dugan,  1996).  The  chemical 
signature  of  landmines  is  affected  by  multiple  environmental  phenomena  that  can  enhance  or  reduce  its  presence  and 
transport,  and  can  affect  the  distribution  of  the  chemical  signature  in  the  environment.  For  example,  the  chemical  can  be 
present  in  the  vapor,  aqueous,  and  solid  phases.  The  distribution  of  the  chemical  among  these  phases,  including  the  spatial 
distribution,  is  key  in  designing  appropriate  detectors,  e.g.  gas,  aqueous  or  solid  phase  sampling  instruments,  and  their 
optimum  use.  A  fundamental  imderstanding  of  the  environmental  conditions  that  affect  the  chemical  signature  is  needed  to 
describe  the  favorable  and  unfavorable  conditions  of  a  chemical  detector  based  survey  to  minimize  the  consequences  of  a 
false  negative.  The  fate  and  transport  of  the  chemical  signature  emanating  fi'om  the  buried  landmine  is  a  fundamental 
property  that  is  poorly  understood.  As  an  initial  step  in  the  evaluation  of  the  landmine  chemical  signature,  a  screening  model 
based  on  pesticide  and  Volatile  Organic  Compound  (VOC)  movement  in  soils  has  been  adapted  to  evaluate  landmine 
chemical  behavior.  This  report  addresses  an  initial  evaluation  of  the  fate  and  transport  of  this  chemical  signature  including 
the  dominant  effects  environmental  conditions  may  have  on  the  success  of  a  chemical  detector  survey.  Future  efforts  to 
develop  more  mechanistic  and  sophisticated  chemical  transport  models  and  the  low  concentration  physical  chemical 
properties  are  needed  to  bridge  the  gap  to  more  realistic  fate  and  transport  conditions. 

Figure  1  shows  a  conceptual  model  of  the  environmental  fate  and  transport  processes  that  impact  the  movement  of  landmine 
chemical  constituents  to  the  surface  for  chemical  detection.  Chemical  vapors  emanate  fi-om  a  buried  landmine  by  permeation 
through  plastic  case  materials  or  through  seals  and  seams,  and  from  the  initial  surface  contamination  of  the  case.  Vapor 
phase  diffusion  transports  molecules  away  from  the  landmine.  The  vapors  may  partition  into  the  aqueous  phase  of  the  soil 
water  which  may  then  be  transported  to  the  surface  through  advection,  driven  by  evapotranspiration  or  to  depth  by 
precipitation  infiltration,  and  through  diffusion  driven  by  concentration  gradients.  Under  extremely  dry  soil  conditions  near 
the  ground  surface,  vapor  phases  may  be  directly  sorbed  to  soil  particles.  When  in  the  liquid  phase,  chemicals  may  also  sorb 
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to  the  soil  particles.  Soil  particle  sorption  can  be  considered  a  temporary  storage  reservoir  for  the  explosive  constituents, 
where  they  may  be  released  under  reversible  partitioning  reactions,  but  some  proportion  may  also  permanently  bound 
through  chemisorption  reactions.  Transformation  and  loss  of  explosive  constituents  also  occurs  during  microbial  degradation 
and  uptake  by  the  roots  of  certain  plant  species. 

2.  CHEMICAL  TRANSPORT  IN  SOILS 


Soils  are  porous  media  with  a  number  of  physico-chemical  properties  that  affect  the  transport  of  explosive  chemicals.  Soil 
bulk  density  is  a  measure  of  the  compaction  of  the  soil  and  is  defmed  as 


[1] 


where  is  the  soil  bulk  density  (g/cm^),  M,  is  the  mass  of  soil  particles  (g),  and  V,  is  the  volume  of  soil  (cm^).  Soils  under 
natural  conditions  have  bulk  densities  ranging  from  1.0  to  1.8  g/cm^.  However,  soils  that  have  been  excavated  and  replaced, 
such  as  during  the  emplacement  of  a  landmine,  may  have  bulk  densities  much  less  than  1 .  The  soil  bulk  density  is  inversely 
proportional  to  the  soil  porosity  as  follows 

PI 


where  is  the  soil  particle  density  (ranges  from  2.6  to  2.8  g/cm^  for  most  soils).  The  soil  porosity,  or  void  volume,  is 
defmed  as 


4>  = 


K+K 

K 


[3] 


where  <j)  is  the  soil  porosity  (emVem^),  is  the  volume  of  soil  water  (cm^)  and  V,  is  the  volume  of  soil  air  (cm^).  Soil 
porosity  values  range  from  0.3  for  sands  to  0.6  for  clay  rich  soils.  The  volumetric  moisture  content  describes  how  much 
water  is  present  in  the  soil  and  changes  greatly  during  precipitation/drainage  events  and  evaporation  conditions.  Volumetric 
water  content  is  defmed  as 


[4] 


where  6  is  the  volumetric  water  content  (cmVcm’).  Soil  moisture  contents  have  values  from  near  zero  up  to  the  soil  porosity 
value.  When  the  soils  are  not  fully  saturated,  the  balance  of  the  soil  pore  space  not  filled  with  water  is  termed  the  air  filled 
porosity,  and  is  defined  as 


[5] 


where  V,  is  the  volumetric  air  content  (cmVcm^).  It  is  often  more  convenient  to  use  soil  saturation  (S,)  because  it  is  a 
measure  of  the  relative  saturation  of  a  particular  soil  pore  space  with  water. 


[6] 


Since  the  explosive  chemicals  can  exist  as  solutes  in  the  soil  water  and  the  movement  of  soil  water  can  be  a  significant 
transport  mechanism,  water  solubility  is  an  important  parameter.  Water  solubility  is  defined  as 


Q  = 


^chem 

V. 


[7] 


where  Cl  is  the  concentration  in  aqueous  phase  (g/cm’  soil  water)  and  M^ion  ts  the  mass  of  chemical  (e.g.  TNT)  (g).  Water 
solubility,  however,  is  not  constant  and  is  typically  an  increasing  function  with  temperature. 


Henry’s  Law  constant  is  a  relative  measure  of  the  amount  of  the  chemical  that  exists  in  the  gas  phase  to  that  in  the  aqueous 
phase  at  equilibrium,  and  is  defined  as 

where  Kh  is  the  Henry’s  Law  constant  (unitless)  and  Cq  is  the  concentration  in  gas  phase  (g/cm3  soil  gas).  Henry’s  Law 
constant  is  also  a  function  of  temperature  because  both  Cq  and  Q  are  functions  of  temperature.  The  soil  partition  coefficient 
is  a  relative  measure  of  how  much  of  the  chemical  is  temporarily  bound  to  the  soil  to  that  in  the  soil  aqueous  phase 

[9] 

where  is  the  linear  soil-water  partition  coefficient  (cm^/g)  for  water  saturated  soils  and  Cj  is  the  concentration  sorbed  on 
the  soil  solid  phase  (g/g  of  soil).  The  soil  water  partition  coefficient  is  often  correlated  with  the  flection  of  organic  carbon 
found  in  the  soils.  In  this  way,  the  variability  between  soils  can  be  reduced.  The  organic  carbon  distribution  coefficient  is 
defined  as 

J  oc 

where  is  the  organic  carbon  distribution  coefficient  and^c  is  the  Section  of  organic  carbon. 


3.  SCREENING  MODELS 

The  environmental  fate  and  transport  of  organic  chemicals  including  volatilization  and  leaching  losses  has  been  used  to 
explore  the  distribution  of  agricultural  pesticides  in  soils  (Mayer  et  al.  1974,  Farmer  et  al.  1980,  and  Jury  et  al.  1980).  These 
models  were  primarily  intended  to  simulate  specific  circumstances.  However,  Jury  et  al.  (1983,  i984a,  1984b,  1984c) 
developed  and  validated  a  general  screening  model  (Behavior  Assessment  Model,  BAM)  that  included  volatilization, 
leaching,  and  degradation  to  explore  the  major  loss  pathways  of  agricultural  pesticides  as  a  function  of  specific 
environmental  conditions.  The  model  simulations  can  be  used  to  assess  the  behavior  of  different  chemicals  under  particular 
environmental  conditions,  but  is  not  intended  to  predict  a  definitive  concentration  distribution  in  the  field.  As  such,  the 
predictions  from  the  screening  model  are  only  an  indication  of  expected  conditions. 


This  model  is  valuable  in  that  it  can  express  the  total  concentration  of  a  chemical  in  the  gas,  aqueous  and  sorbed  phases.  The 
total  concentration  is  expressed  as 


”  Pb^S 

where  Cj  is  the  concentration  sorbed  to  the  soil,  Q  is  the  solute  concentration  in  the  aqueous  phase,  and  Cq  is  the  gas  phase 
concentration.  In  addition,  Jury  (1983)  shows  how  equation  [1 1]  can  be  rewritten  in  terms  of  one  of  the  variables  alone 


where 


Cj.  =  RsCs  =  RlCj^  =  RgCc 

[12] 

e  K„ 

=  Pb  Kj 

[13] 

=  Pi,Kj  +6  +aK„,  and 

[14] 

^G-Pb  rr  + 

[15] 

are  the  solid,  liquid  and  gas  phase  partition  coefficients,  respectively. 

An  adaptation  of  the  BAM  w'as  developed  to  be  applicable  to  the  conditions  of  contaminated  soil  buried  under  a  known 
depth  of  clean  soil  -  Buried  Chemical  Model,  BCM  (Jury  et  al.,  1990).  Simulations  based  on  a  modification  of  Jury’s  BCM 
are  used  in  this  report  to  simulate  the  behavior  of  the  chemical  signature  from  buried  landmines.  The  Buried  Chemical 
Model  of  Jury  et  al.  (1990)  is  based  on  the  following  assumptions.  A  detailed  discussion  of  these  assumptions  is  given  in 
Jury  et  al.  (1990). 

1.  The  chemical  may  adsorb  on  the  solid  phase,  be  dissolved  in  the  aqueous  phase,  or  exist  in  the  vapor  phase. 

2.  The  chemical  flux  is  the  sum  of  the  vapor  flux  and  the  dissolved  solute  flux. 

3.  The  porous  medium  factors  for  gas  and  liquid  phase  diffusion  are  given  by  the  Millington  and  Quirk  (1961)  model  as 

extended  for  liquid  diffusion  by  Jury  et  al.  (1983). 

4.  The  chemical  will  undergo  first-order  degradation  due  to  biological  and  chemical  effects, 

5.  Chemical  movement  is  one  dimensional. 

6.  The  adsorbed  and  dissolved  phases  undergo  reversible,  linear  adsorption. 

7.  The  dissolved  and  gaseous  phase  concentrations  are  related  through  Henry’s  law. 

8.  The  soil  properties  are  constant  in  space  and  time. 

9.  Water  flux  is  constant  in  space  and  time  (relaxed  in  the  present  application). 

10.  Volatilization  of  the  chemical  to  the  atmosphere  is  by  vapor  diffusion  through  an  air  boundary  layer  of  constant 
thickness. 


In  the  present  implementation  of  Jury’s  model,  a  constant  source  term  has  been  added  to  reflect  the  chemical  source  from  the 
landmine  at  a  specific  location. 


Under  these  assumptions  (including  the  source  term)  the  model  formulation  becomes 


dCj 

dt 


+  \xCj  — 


d^Cr 


dz^  ^  dz 


[16] 


where  Cr  is  the  total  chemical  concentration,  p  is  the  biochemical  decay  constant,  and  a  is  the  source  term.  The  effective 
velocity  (V^)  is  defined  as 


[17] 


F.= 


J. 


Ph^d  +6 


where  J,  is  the  precipitation/evaporation  flux.  The  effective  diffusion  coefficient  (D^)  of  the  chemical  is  defined  as 

r\o  I 
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[18] 


where  Z)“  is  the  diffusivity  of  the  gas  phase  of  the  chemical  in  air  and  D"  is  the  diffusivity  of  the  chemical  in  aqueous 

phase .  The  boundaiy  conditions  for  the  problem  are  diffusion  through  a  boundary  layer  at  the  upper  surface,  and  a  zero 
chemical  concentration  at  infinity  at  the  lower  boundary.  These  boundary  conditions  can  be  expressed  as 


-A 


da 


dz 
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[20] 


[21] 


=  0 

The  initial  conditions  are  an  initial  concentration,  Co,  over  an  interval  from  L  to  W,  or 

Q(2,0)  =  0  L)z)W 

Q(z,0)  =  Co  L<z<W 

The  above  model  without  the  source  term  results  in  a  closed  form  solution  as  a  function  of  space  and  time;  the  results  are 
rather  lengthy  and  will  not  be  presented  here  but  are  given  by  Jury  et  al.  (1983, 1990).  In  the  present  simulations,  the 
assumption  of  constant  water  flux  in  time  is  relaxed.  Therefore,  sequences  of  water  fluxes  representing  desired  conditions 
(rainfall  followed  by  evaporation)  can  be  simulated  to  determine  the  effect  of  water  flux  variations  on  the  location  of  TNT  in 
the  soil  and  the  surface  TNT  vapor  flux.  A  numerical  solution  was  developed  and  verified  by  comparison  to  the  results  given 
by  Spencer  et  al.  (1988)  and  Jury  et  al.  (1990)  (Phelan  and  Webb,  1997). 


4.  LANDMINE  SOURCE  TERM 

The  total  mass  of  the  initial  deposit  of  chemicals  after  first  emplacement  of  the  landmine  in  soils  is  critical  in  the  estimation 
of  soil  concentrations  and  surface  vapor  fluxes  (Phelan  and  Webb,  1988).  For  the  model  used  in  this  analysis,  it  has  been 
assmned  that  the  entire  surface  contamination  was  completely  and  uniformly  transferred  to  the  soil  just  prior  to  the  beginning 
of  the  simulation  runs.  Surface  contamination  data  (Hogan  et  al.,  1992)  showed  a  median  surface  contamination  of  15 
ng/cm^  from  42  domestic  and  foreign  landmines.  However,  the  surface  contamination  values  ranged  from  below  instrument 
detection  limits  to  300-500  ng/cm2  and  some  outliers  as  high  as  1000-5000  ng/cm2.  Using  the  median  surface 
contamination  and  the  dimensions  of  an  anti-tank  (AT)  mine  of  30  cm  diameter  by  10  cm  high,  the  surface  contamination 
would  provide  3.5x10'’  g  of  TNT  for  initial  distribution  in  the  soil.  Using  the  volume  of  the  AT  mine  that  this  mass  of  TNT 
is  distributed  into,  the  initial  concentration  (CJ  would  be  5x10’’  pg/cm’. 

The  continuous  release  of  chemicals  by  permeation  or  leaks  through  seals  and  seams  is  also  important  if  the  rate  is 
significantly  large  enough.  However,  few  data  exist.  One  effort  established  constant  source  term  emanation  rates  derived 
from  vapor  collection  chamber  experiments  on  two  mines  (Spangler,  1975).  Values  ranged  from  10''*  to  10''*  g/cm^-s.  The 
higher  rate  of  10'”  g/cm^-s  (8.6x10'*  pg/cm^-day)  was  used  in  these  simulations. 


5.  PHYSICAL  CHEMICAL  PROPERTIES  OF  EXPLOSIVES 


The  principal  explosive  chemicals  found  in  landmines  are  TNT  and  RDX  (NGIC,  1995).  DNT,  as  a  production  by-product 
of  TNT,  is  also  considered  to  be  a  significant  signature  chemical  for  buried  landmines.  As  a  group,  these  chemicals  have 
very  low  vapor  densities  and  moderately  low  water  solubility’s.  Table  1  shows  these  properties  and  the  Henry’s  Law 
constant  at  20‘’C  (Phelan  and  Webb,  1997). 


Table  1.  Vapor  Density,  Water  Solubility  and  Henry’s  Law  Constant  of  Explosive  Compounds  at  20°C 

_ TNT  DNT  RDX 

Vapor  Density  (pg/m^  43.5  122  0.024 

Water  Solubility  (mg/1)  130  270  50 

Hemy’s  Law  Constant _ 3.35E-7  4.51E-7  4.73E-10 


Pennington  and  Patrick  (1990)  measured  the  soil  water  partition  coefficient  (K,,)  of  TNT  in  fourteen  soils  from  locations 
across  North  America.  The  mean  value  was  3.8  cmVg  with  a  standard  deviation  of  1.34.  The  highest  value  was  6.8  cmVg 
and  the  lowest  value  was  2.3  cmVg.  Xue  et  al.  (1995)  evaluated  two  soils  and  showed  mean  values  for  TNT  of  2.66  cmVg 
and  3.64  cmVg.  DNT  and  RDX  have  very  little  data.  Xiie  et  al.  (1995)  showed  values  for  RDX  of  1.59  cmVg  and  1.57 
cmVg.  McGrath  (1995)  showed  a  value  of  251  for  DNT.  For  the  fourteen  soils  evaluated  by  Pennington  and  Patrick 
(1990),  the  mean  value  for  the  fraction  of  organic  carbon  was  0.0173  with  a  standard  deviation  of  0.01 1.  Using  these  values, 
the  Kj  for  DNT  has  a  mean  value  of  4.4±2.7  cmVg  (one  std.  dev.).  In  summary,  the  soil  water  partitioning  coefficients  for 
TNT,  DNT  and  RDX  all  fall  into  an  approximate  range  between  1.5  and  7.0  cmVg.  This  is  a  rather  narrow  range  as  common 
chemicals  can  have  values  one  to  two  orders  of  magnitude  lesser  and  greater  than  these. 

The  biochemical  half-life  of  explosives  in  near  surface  soils  has  not  been  studied  well  outside  of  the  biotreatment  technology 
area  for  contaminated  soils.  However,  long-term  surface  soil  degradation  tests  at  Los  Alamos  National  Laboratory  (LANL) 
followed  the  degradation  of  soils  doped  with  1000  mg/kg  of  various  explosives  over  20  years  (Dubois  and  Baytos,  1991). 
Table  2  shows  the  half-lives  estimated  from  these  long-term  experiments. 

Table  2.  Estimates  of  Half-Lives  of  Explosives  from  LANL  Long-Term  Surface  Soil  Tests 


Explosive 

Half-Life  (years) 

TNT 

1 

RDX 

36 

HMX 

39 

PETN 

92 

6.  SIMULATION  PARAMETERS 


Using  the  modified  Buried  Chemical  Model  (BCM),  simulations  were  performed  using  a  landmine  that  has  contributed  an 
initial  soil  concentration  (C^)  based  on  the  median  surface  contamination  of  the  landmine  (Hogan  et  al.,  1992)  and  the 

constant  flux  (JJ  from  Spangler,  1975.  The  diffiisivity  of  gas  in  air  (  Z)” )  and  diffusivity  of  liquid  in  water  (  Df  )  were 

selected  from  Jury  et  al.  (1983).  The  biochemical  half-life  value  of  365  days  was  selected  from  a  long  term  field  experiment 
(Dubois  and  Bayton,  1991).  The  physico-chemical  properties  of  TNT  were  used  at  25®C. 

The  precipitation/evaporation  rates  and  periods  followed  in  several  of  the  simulations  here  were  the  low  desert  scenario  from 
Phelan  and  Webb  (1997).  This  scenario  was  derived  from  data  found  in  HELP  (Hydrological  Evaluation  of  Landfill 
Performance)  model  (Schroeder  et  al.,  1994a  and  1994b).  The  HELP  model  showed  that  the  low  desert  had  1  day  of 
precipitation  followed  by  7  days  of  evaporation.  For  simplicity,  total  precipitation  and  total  evaporation  for  each  cycle  are 
assumed  to  be  equal  and  for  these  simulations  the  cycles  were  continued  for  approximately  four  years.  Table  3  shows  the 
input  parameters  used  in  the  simulations. 


Table  3.  Simulation  Parameters 


parameter 

units 

base  case 

variant  cases 

e 

cmV  cm^ 

0.20 

Figure  2:  0.243,0.287 

<!> 

cmV  cm^ 

0.434 

* 

pb 

g/cm^ 

1.5 

Figure  2;  1.25,1.0 

Kd 

cmVg 

3.8 

* 

Kh 

5.9E-7 

air  boundary 
layer 

cm 

0.5 

♦ 

^1/2 

days 

365 

♦ 

C„ 

Hg/cm’ 

5E-3 

*  “ 

Jc 

pg/cm^-day 

8.6E-6 

* 

d; 

cmVday 

0.432 

♦ 

cmVday 

4320 

* 

burial  depth, 
top 

cm 

10 

Figure  3:  5, 1,  0 

Figure  6,  7,  8  &  9:  5 

burial  depth, 
bottom 

cm 

20 

Figure  3:  15,11,10 

Figure  6,  7,  8&9:  15 

precipitation 

days 

1 

varies 

evaporation 

days 

7 

varies 

precipitation 

rate 

cm/day 

0.44 

Figure  4,  5:  0 

evaporation 

rate 

cm/day 

-  0.063 

Figure  4,  5:  0 

*  -  same  as  the  base  case 


7.  DISCUSSION 

Initial  application  of  the  screening  model  (Phelan  and  Webb,  1997)  showed  that  explosive  compounds  such  as  TNT,  DNT 
and  RDX  will  have  over  90%  of  the  mass  fraction  sorbed  to  the  soil  solid  phase,  up  to  1 0%  present  in  the  soil  aqueous  phase 
and  less  than  lE-6%  in  the  soil  vapor  phase.  The  implications  for  this  are  that  transport  of  these  chemicals  in  soils  will  be 
dominated  by  movement  as  a  solute  in  the  aqueous  phase.  In  addition,  at  soil  water  contents  above  about  0.1  cmVcm^,  the 
effective  diffiisivity  is  dominated  by  liquid  diffusion  by  one  to  five  orders  of  magnitude.  Next,  the  screening  model  was 
applied  to  three  different  environmental  scenarios  and  results  showed  that  maximum  surface  vapor  fluxes  were  very  low, 
requiring  a  sampling  and  concentration  factor  of  10^  using  currently  available  laboratory  instrument  detection  limits.  Surface 
soil  concentrations  showed  that  sampling  modest  amounts  of  soil  would  require  a  preconcentration  factor  of  about  2. 

Next,  the  impact  of  variations  in  several  of  the  various  input  parameters  was  explored  (Phelan  and  Webb,  1998).  It  was 
found  that  the  variations  in  the  Henry’s  Law  constant  are  directly  proportional  to  the  changes  seen  in  the  steady-state  surface 
vapor  flux.  Small  changes  in  the  soil  water  partitioning  coefficient  made  large  changes  in  the  lag  period  and  steady-state 
surface  vapor  flux.  The  initial  source  term  and  continuous  source  flux  of  chemicals  from  the  landmines  is  poorly  imderstood 
and  expected  to  be  highly  variable.  These  simulations  found  that  with  the  baseline  case,  if  the  continuous  source  flux  was 
absent,  there  was  no  significant  difference  in  the  surface  vapor  flux  or  subsurface  distribution  at  the  end  of  the  simulation 
period.  This  implies  that  the  continuous  source  flux  may  be  much  less  important  than  the  initial  surface  contamination.  It 
appears  that  the  magnitude  of  the  surface  vapor  flux  is  directly  proportional  to  the  amount  of  the  initial  surface 
contamination.  The  biochemical  half-life  is  another  parameter  that  is  likely  to  have  veiy  different  values  depending  on  the 
location  and  climatic  conditions.  Simulations  over  ten  years  showed  that  the  steady  state  surface  vapor  flux  declines  steadily 
when  the  biochemical  half-life  becomes  smaller  than  one  year.  Finally,  the  importance  of  heavy  precipitation  (such  as  a 
monsoon  season)  followed  by  a  dry  season  was  explored.  The  impact  of  the  heavy  precipitation  was  to  lower  the  surface 


flux  seven  orders  of  magnitude;  however,  the  evaporation  period  that  followed  returned  the  surface  flux  to  approximately 
the  pre-monsoon  surface  flux. 


This  current  effort  has  evaluated  the  impact  of  changes  to  the  soil  bulk  density,  burial  depth,  concentration  in  the  soil  layer, 
water  flux,  and  the  conditions  that  promote  an  enhanced  surface  layer.  When  a  landmine  is  placed,  the  soil  surrounding  the 
mine  is  not  compacted  to  the  original  bulk  density.  Over  time  the  soil  bulk  density  will  increase.  Figure  2  shows  that  soils 
with  smaller  soil  bulk  densities  will  have  a  shorter  lag  period  and  a  steady  state  surface  flux  about  two  orders  of  magnitude 
greater  than  the  base  case. 
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Figure  2.  Effect  of  Soil  Bulk  Density  Changes  on  Surface  Flux 

The  effect  of  burial  depth  appears  to  be  a  very  critical  parameter,  even  over  small  distances  of  a  few  centimeters.  Figure  3 
shows  how  the  lag  time  for  the  surface  vapor  flux  becomes  dramatically  shorter  by  moving  the  top  of  the  initial  source  zone 
up  from  10  cm  to  5  cm  to  1  cm  and  0  cm. 

1.00E*«J 
1.00e4C 
i.ooe-04 

t.SOE-06 

vooe-w 

t.OOE-tO 

tjOOE<12 
^  t.(IOE-1<< 

•  i.ooE-ie 

i.fiOE-ie 

i.ooe*2D 

1.00E-22 

1.00E-24 

Figure  3.  Effect  of  Burial  Depth  on  Surface  Flux 

An  important  assumption  made  in  the  current  numerical  application  to  the  landmine  problem  is  the  assumption  of  a  uniform 
explosive  concentration  between  the  top  and  bottom  of  the  mine.  In  reality,  the  explosive  is  concentrated  on  the  top  and  the 
bottom  surfaces,  the  mine  is  an  impediment  to  transport,  and  the  mine  is  not  porous  media.  In  order  to  assess  the  implication 
of  smearing  the  concentration  over  the  entire  depth,  Jury’s  original  BCM  without  a  source  term  and  with  zero  water  flux  has 
been  applied  to  the  case  of  two  buried  chemical  sources.  The  top  source  started  at  the  top  of  the  landmine  location,  while  the 
bottom  source  ended  at  the  bottom  of  the  landmine,  with  a  zero  concentration  region  in  the  center  of  the  landmine  location. 
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A  constant  mass  was  used,  and  the  concentration  changed  with  the  source  width,  which  varied  from  5.0  cm,  the  base  case, 
down  to  0.0005  cm.  The  5.0  cm  base  case  results  in  a  uniform  initial  concentration  between  the  top  and  bottom  of  the  mine, 
while  the  0.0005  cm  case  approximates  two  sources  at  the  top  and  bottom  surfaces.  Fig  4  shows  the  vapor  flux  at  the  ground 
surface  as  a  function  of  time  for  the  various  source  widths.  The  parameters  in  this  case  are  the  same  as  for  the  base  case  in 
Table  3  without  a  source  term  and  with  no  precipitation/evaporation.  The  results  are  similar,  though  not  identical,  to  results 
given  in  Figure  2  (rho  b  =  1.5). 
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Figure  4.  Effect  of  Layer  Concentration  on  Surface  Vapor  Flux 

Even  if  the  parameters  were  the  same,  some  differences  would  exist,  particularly  at  earlier  times  and  at  lower  concentrations, 
due  to  numerical  issues  such  as  roundoff  and  numerical  diffusion.  From  Figure  4,  the  vapor  flux  increases  as  the  width 
decreases,  which  is  expected  because  a  reduction  in  the  width  essentially  moves  the  top  source  closer  to  the  ground  surface. 
Figure  5  shows  the  ratio  of  the  vapor  flux  from  the  0.0005  cm  width  case  to  the  5  cm  case.  The  ratio  decreases  with  time 
from  about  50  to  a  value  of  10  at  late  times. 
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Figure  5.  Ratio  of  Surface  Flux  from  a  0.0005  cm  Layer  to  5  cm  Layer 

In  addition  to  the  effect  of  soil  moisture  content,  the  effect  of  precipitation  (positive  water  flux,  JJ  and  evaporation  (negative 
water  flux,  JJ  is  probably  one  of  the  most  important  environmental  factors  in  the  transport  of  explosive  chemicals  in  soils. 
Figure  6  shows  that  with  only  precipitation  occurring,  the  surface  flux  is  about  3  orders  of  magnitude  less  than  the  case  of 
zero  precipitation  or  evaporation.  The  case  of  constant  evaporation  is  about  2  orders  of  magnitude  greater  than  the  zero 
water  flux  case.  If  one  examines  the  model  formulation,  the  mass  transport  upwards  is  controlled  by  the  effective  diffusion 
(Dg)  and  the  effective  chemical  velocity  (Vg).  In  the  constant  precipitation  case,  upward  mass  transport  is  a  function  of  Dj 


minus  V^.  In  the  zero  precipitation/evaporation,  upward  mass  transport  is  a  function  of  only  Dg.  In  the  constant  evaporation 
case,  upward  mass  transport  is  a  function  of  Dg  plus  V^. 


Figure  6.  Effect  of  Water  Flux  (Precipitation/Evaporation)  on  Surface  Flux 

The  occurrence  of  a  surface  soil  layer  that  is  greater  in  concentration  than  the  subsurface  soil  layers,  or  a  surface  “crust”,  has 
been  observed  and  modeled  with  agricultural  pesticides  (Spencer  et  al.,  1988).  This  type  of  behavior  is  also  thought  to  occur 
with  explosive  compounds  due  to  similarities  in  the  physical/chemical  properties  and  some  evidence  from  field  surveys  and 
lab  experiments.  Simulation  runs  were  completed  to  evaluate  what  influences  the  creation  of  the  enhanced  concentrations  in 
the  surface  soil  layers.  Initial  simulation  runs  (Phelan  and  Webb,  1997)  used  cyclic  precipitation/evaporation  that  was  equal 
in  magnitude.  This  condition  did  not  create  an  enhanced  surface  layer.  In  order  to  create  an  enhanced  surface  layer,  enough 
of  the  mass  must  be  transported  from  deeper  regions  to  the  ground  surface.  This  condition  only  occurs  during  evaporation 
conditions  and  in  Figure  7  the  buried  chemical  layer  is  shown  to  move  upward  until  it  intersects  with  the  ground  surface. 
Figure  8  shows  more  detail  on  the  depth  and  magnitude  of  the  enhanced  layer.  It  is  believed  that  the  air  boundary  layer  and 
the  low  Henry’s  Law  Constant  (Kh)  contribute  to  the  formation  of  the  enhanced  surface  layer  (Spencer  et  al.,  1988).  It 
appears  that  the  upward  transport  through  the  soil  exceeds  the  loss  through  the  air  boundary  layer.  Transport  though  the  air 
boundary  layer  is  controlled  by  vapor  diffusion  and  limited  by  the  transfer  of  chemical  from  the  aqueous  phase  to  the  gas 
phase  by  the  very  low  Kh- 
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Figure  7.  Upward  Transport  and  Development  of  a  Surface  Layer 

Figure  9  shows  the  effect  of  developing  the  enhanced  surface  layer  with  60  days  of  evaporation  (-0.5  cm/day),  followed  by 
precipitation  for  5  days  (0.5  cm/day).  The  enhanced  surface  layer  found  in  the  top  0.1  cm  of  soil  is  transported  down  leaving 
just  a  small  enhancement  at  a  depth  of  about  0.5  cm.  Another  simulation  was  run  that  included  the  same  evaporation  and 


precipitation,  but  was  followed  by  another  5  day  evaporation  period  (-0.5  cm/day)  and  the  surface  enhancement  returned  at 
about  the  same  concentration. 


Figure  8.  Detail  of  the  Surface  Layer  Formed  by  Evaporation 


Figure  9.  Effect  of  Precipitation  on  the  Enhanced  Surface  Layer 
8,  SUMMARY  AND  CONCLUSIONS 

Two  previous  efforts  have  evaluated  the  environmental  fate  and  transport  of  chemical  signature  molecules  from  buried 
landmines  (Phelan  and  Webb,  1997  and  1998).  This  effort  evaluated  the  effect  of  changes  in  soil  bulk  density,  burial  depth, 
source  term  thickness,  precipitation  or  evaporation,  and  the  elements  that  produce  an  enhanced  surface  layer.  Smaller  soil 
bulk  densities  change  the  soil  porosity  and  moisture  content  at  constant  soil  saturations.  Lower  soil  bulk  densities  increase 
the  rate  of  transport  as  indicated  by  a  shorter  lag  time  and  reach  a  higher  steady  state  surface  vapor  flux.  Burial  depth  has  a 
dramatic  effect  on  the  lag  time.  With  just  a  one  centimeter  soil  cover,  there  appears  to  be  almost  no  lag  time  and  the  surface 
vapor  flux  is  five  orders  of  magnitude  greater  than  a  soil  cover  of  ten  centimeters.  The  assumption  of  a  uniform 
concentration  over  the  depth  of  the  landmine  leads  to  under  prediction  of  the  vapor  flux  at  the  land  surface  by  about  an  order 
of  magnitude  or  more.  The  appearance  of  an  enhanced  surface  layer  is  the  result  of  evaporation  driven  mass  transport  that  is 
constrained  by  vapor  diffusion  through  the  air  boundary  at  the  ground  surface  and  the  low  aqueous  to  vapor  transfer  as  a 
function  of  the  low  Henry’s  Law  Constant.  However,  there  are  probably  many  other  uncertainties  in  the  present  model 
which  affect  the  quantitative  results  as  well.  The  one-dimensional  assumption  demands  that  the  mine  itself  is  permeable,  and 
the  assumption  of  constant  liquid  content,  both  spatially  and  temporally,  is  obviously  a  great  simplification.  In  order  to 
address  these  and  other  issues,  a  multidimensional  mechanistic  code  is  being  developed  for  application  to  the  landmine 


problem.  This  code,  which  is  based  on  the  TOUGH  code  from  Lawrence  Berkeley  Laboratory  (Pruess,  1987,1991),  will 
consider  air,  water  vapor,  and  explosive  vapor  mass  and  heat  flow  in  a  porous  media  and  will  be  able  to  address  many  of 
these  questions. 
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ABSTRACT 

The  detection,  and  removal,  of  buried  landmines  is 
one  of  the  most  important  international  problems 
facing  the  world  today.  Numerous  detection 
strategies  are  being  developed,  including  infrared, 
electrical  conductivity,  ground«penetrating  radar,  and 
chemical  sensors.  Chemical  sensors  rely  on  the 
detection  of  TNT  molecules  from  buried  landmines, 
which  are  transported  by  advection  and  diffusion  in 
the  soil.  As  part  of  this  effort,  numerical  models  are 
being  developed  to  predict  TNT  transport  in  soils 
including  the  effect  of  precipitation  and  evaporation. 
Modifications  will  be  made  to  TOUGH2  for 
application  to  the  TNT  chemical  sensing  problem. 
The  understanding  of  the  environmental  fate  and 
transport  of  TNT  in  the  soil  will  affect  the  design, 
performance  and  operation  of  chemical  sensors  by 
indicating  preferred  sensing  strategies, 

INTRODUCTION 

The  goal  of  locating  buried  landmines  is  a 
significant  challenge  to  science  and  technology.  The 
chemical  signature  of  landmines  is  affected  by 
multiple  environmental  phenomena  that  can  enhance 
or  reduce  its  presence  and  transport  behavior,  and  can 
affect  the  distribution  of  the  chemical  signature  in  the 
environment.  For  example,  the  chemical  can  be 
present  in  the  vapor,  aqueous,  and  solid  phases.  The 


distribution  of  the  chemical  among  these  phases, 
including  the  spatial  distribution,  is  key  in  designing 
appropriate  detectors,  e.g.,  gas,  aqueous  or  solid 
phase  sampling  instruments.  A  fundamental 
understanding  of  the  environmental  conditions  that 
affect  the  chemical  signature  is  needed  to  describe 
the  favorable  and  unfavorable  conditions  of  a 
chemical  detector  based  survey  to  minimize  the 
consequences  of  a  false  negative. 

The  fate  and  transport  of  the  chemical  signature 
emanating  from  the  buried  landmine  is  poorly 
understood.  As  an  initial  step  in  the  evaluation  of  the 
landmine  chemical  signature,  a  screening  model 
based  on  pesticide  and  Volatile  Organic  Compound 
(VOC)  movement  in  soils  has  been  adapted  to 
evaluate  landmine  chemical  behavior.  Future  efforts 
to  develop  more  mechanistic  and  sophisticated 
chemical  transport  models  are  needed  to  bridge  the 
gap  to  more  realistic  fate  and  transport  conditions. 

Figure  1  shows  a  conceptual  model  of  the 
environmental  fate  and  transport  processes  that 
impact  the  movement  of  landmine  chemical 
constituents,  such  as  TNT  and  DNT,  to  the  land 
surface  for  chemical  detection.  Chemical  vapors 
emanate  from  a  buried  landmine  by  permeation 
through  plastic  case  materials  or  ‘"leakage”  through 
seals  and  seams,  and  from  surface  contamination  of 
the  case. 


Figure  1,  Environmental  Fate  and  Transport  Model  for  Chemical  Detection  of  Buried  Landmines. 
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Table  1.  TNT  and  DNT 

Properties  at 

20''C 

TNT 

DNT 

Vapor  Density  (pg/m'^) 

43.5 

122 

Water  Solubility  (mg/l) 

1.30 

270 

Henry’s  Law  Constant 

3.35E-7 

4.51E-7 

Sorption  Coefficient  (cm*Vg) 

3.8 

4.4 

CHEMICAL  PROPERTIES 
The  chemical  properties  of  TNT  and  DNT  are 
important  in  determining  the  transport  rate  of  these 
vapors  through  the  soil.  These  chemical  vapors  exist 
in  the  gas,  liquid,  and  solid  phases  of  the  soil. 
Typical  properties  for  TNT  and  DNT  are  shown  in 
Table  1.  Because  of  the  low  value  of  Henry’s 
constant  and  the  value  of  the  soil  water  partition 
coefficient,  about  90%  of  the  explosive  mass  fraction 
is  sorbed  to  the  soil  solid  phase,  about  10%  is  in  the 
water;  and  less  than  10’^%  is  in  the  gas  phase  as 
shown  in  Figure  2. 
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Figure  2.  Phase  Mass  Fraction  of  TNT 

The  transport  rate  of  TNT  in  soil  can  be  estimated 
by  evaluating  pure  diffusion  conditions.  An  effective 
diffusivity  can  be  defined  for  the  total  chemical 
concentration  by  considering  the  distribution  among 
the  phases  (Jury  et  al.,  1983).  By  applying  the 
Millington  and  Quirk  (1961)  tortuosity  relationship  to 
the  liquid  phase,  the  effective  diffusivity  for  the  total 
chemical  concentration  can  be  expressed  as 


(1) 


Figure  3  shows  the  variation  in  this  effective 
diffusivity  with  water  content.  Note  that  the 
diffusivity  value  is  low  due  to  the  value  of  Henry’s 
constant  and  sorption  onto  the  solid  phase  that  acts  as 


a  sink  for  the  explosive  chemical.  The  effective 
diffusivity  is  generally  much  higher  al  higher 
moisture  contents. 


Figure  3,  TNT  Effective  Diffusivity 


INITIAL  TRANSPORT  STUDIES 

Jury  and  his  colleagues  (Jury  et  al.,  1983, 
19841,b,c)  developed  a  one-dimensional  screening 
model  to  study  the  behavior  of  various  pesticides 
under  different  environmental  conditions.  This 
model  was  subsequently  extended  to  buried 
chemicals,  such  as  VOCs,  by  Jury  et  al.  (1990).  TNT 
properties  (Henry’s  constant  and  sorption  coefficient) 
are  very  similar  to  some  pesticides,  especially 
Prometron.  Therefore,  Jury’s  screening  model  has 
been  used  for  some  initial  studies  of  the  TNT 
transport  in  soils  from  buried  landmines. 

In  addition,  Prometron  exhibits  some  interesting 
behavior  that  may  be  particularly  important  for  the 
sensing  of  TNT  for  landmine  detection.  Under 
evaporation  conditions,  a  surface  “crust”,  or  a  soil 
layer  that  is  greater  in  concentration  than  the 
subsurface  soil,  has  been  observed  in  laboratory  tests; 
this  surface  “crust”  is  also  predicted  by  Jury’s 
screening  model  (Spencer  et  al.,  1988).  Some 
evidence  of  this  type  of  behavior  for  TNT  has  also 
been  noted  in  some  field  surveys  and  lab 
experiments,  although  the  data  are  not  definitive. 
The  occurrence  of  a  surface  “crust”  would  greatly 
enhance  the  concentration  available  to  chemical 
sensors  and  the  efficiency  of  the  technique. 

Initial  studies  of  the  transport  rate  of  TNT  in  soils 
from  landmines  have  been  conducted  by  Phelan  and 
Webb  (1997,  1998a,b)  using  Jury’s  model.  The 
results  of  Phelan  and  Webb  for  landmine  detection 
indicated  a  significant  influence  of  the  soil  type  and 
environmental  conditions,  including  precipitation  and 
evaporation,  on  the  TNT  flux  at  the  soil  surface 
which  is  available  to  chemical  detectors. 


Surface  Flux  (ug/cm2-day) 


Table  2.  Phase  Specific  Concentration  of  TNT  at  the  Ground  Surface  After  One  Year 


Volumetric  Water  Content/Saturation 


Concentration 

Units 

0.20/0.46 

0.25/0.58 

0.30/0.69 

Solid  Phase 

fig  TNT/g  soil 

1.8E-8 

3.1E-6 

2.8E-5 

Liquid  Phase 

fig  TNT/ml  soil  water 

4.8E-9 

8.4E-7 

7.6E-6 

Gas  Phase 

fig  TNT/cm^  soil  air 

2.8E-15 

5.0E-13 

4.5E-12 

Results  from  this  screening  model  are  shown  below 
for  a  landmine  buried  from  5  cm  to  15  cm  beneath 
the  surface;  details  are  given  by  Phelan  and  Webb 
(1997,  1998a,b).  Note  that  the  screening  model  was 
developed  to  assess  the  behavior  of  different 
chemicals  under  specific  environmental  conditions;  it 
is  not  intended  as  a  purely  predictive  model  due  to  a 
number  of  simplifying  assumptions,  such  as  constant 
soil  moisture  content.  Therefore,  these  results  are 
only  an  indication  of  expected  conditions,  and  more 
detailed  numerical  methods,  such  as  TOUGH2,  are 
necessary  for  a  fully  predictive  simulation. 

Figure  4  shows  the  TNT  surface  flux  at  the  land 
surface  for  a  Gulf  coastal  lowlands  soil  type  as  a 
function  of  soil  saturation.  The  oscillations  in  the 
surface  flux  are  a  result  of  precipitation/evaporation 
cycles,  which  were  constant  over  the  simulation.  As 
the  soil  saturation  increases,  the  surface  flux 
increases  dramatically. 


Time  (days) 

Figure  4.  Surface  Flux  of  TNT 

Figure  5  shows  the  surface  distribution  of  TNT 
after  1  year  for  the  three  soil  saturations.  The 
landmine  was  buried  from  5  to  15  cm  below  the 
ground  with  an  initial  concentration  based  on 
contamination  on  the  landmine  casing.  In  addition, 
landmines  “leak”  TNT  through  the  casing,  which  was 
represented  by  a  source  at  10  cm.  The  movement  of 
the  TNT  away  from  the  landmine  is  slow  and  is  a 
function  of  the  liquid  saturation. 


Note  the  low  total  concentrations  in  Figure  5. 
These  total  concentrations  can  be  further  broken 
down  into  solid,  liquid,  and  gas  phase  values  as 
summarized  in  Table  2,  which  indicate  the  extremely 
small  concentrations  available  in  the  gas  phase.  This 
information  will  be  valuable  in  the  design  and 
operation  of  chemical  sensors  for  landmine  detection. 


Figure  5.  Subsurface  Distribution  of  TNT 
After  One  Year 

The  effect  of  burial  depth  is  a  critical  parameter. 
Figure  6  shows  how  the  lag  time  for  the  surface 
vapor  flux  becomes  dramatically  shorter  by  moving 
the  top  of  the  initial  source  zone  up  from  1 0  cm  to  5 
cm  to  1  cm  and  0  cm. 


Time  (days) 

Figure  6.  Effect  of  Burial  Depth  on  Surface  Flux 


Surface  Flux  (ug/cm2»day) 


The  effect  of  precipitation  (positive  water  flux,  Jw) 
and  evaporation  (negative  water  flux,  is  one  of  the 
most  important  environmental  factors  in  the  transport 
of  explosive  chemicals  in  soils.  Figure  7  shows  that 
with  only  precipitation  occurring,  the  surface  flux  is 
about  3  orders  of  magnitude  less  than  the  case  of  zero 
precipitation  or  evaporation.  The  case  of  constant 
evaporation  is  about  2  orders  of  magnitude  greater 
than  the  zero  water  flux  case.  If  one  examines  the 
model  formulation,  the  mass  transport  upwards  is 
controlled  by  the  effective  diffusion  (De)  and  the 
effective  chemical  velocity  (Ve).  In  the  constant 
precipitation  case,  upward  mass  transport  is  a 
function  of  De  minus  Ve.  For  the  case  of  zero 
precipitation/evaporation,  upward  mass  transport  is  a 
function  of  only  De-  In  the  constant  evaporation  case, 
upward  mass  transport  is  a  function  of  De  plus  Ve. 


Figure  7.  Effect  of  Water  Flux 
(Precipitation/Evaporation)  on  Surface  Flux 

The  occurrence  of  a  surface  soil  layer  that  is  greater 
in  concentration  than  the  subsurface  soil  layers,  or  a 
surface  “crust”,  was  discussed  earlier.  Simulations 
were  performed  to  evaluate  what  influences  the 
creation  of  the  enhanced  concentrations  in  the  surface 
soil  layers.  Initial  simulation  runs  (Phelan  and  Webb, 
1997,  1998a)  used  cyclic  precipitation/evaporation 
that  was  equal  in  magnitude;  this  condition  did  not 
create  an  enhanced  surface  layer.  In  order  to  create 
an  enhanced  surface  layer,  enough  of  the  mass  must 
be  transported  from  deeper  regions  to  the  ground 
surface.  This  condition  only  occurs  during 
evaporation  conditions.  In  Figure  8  the  buried 
chemical  layer  is  shown  to  move  upward  until  it 
intersects  with  the  ground  surface.  Figure  9  shows 
the  depth  and  magnitude  of  the  enhanced  layer.  It  is 
believed  that  the  air  boundary  layer  and  the  low 
Henry’s  Law  Constant  (Kh)  contribute  to  the 
formation  of  the  enhanced  surface  layer  (Spencer  et 


al.,  1988).  It  appears  that  the  upward  transport 
through  the  soil  exceeds  the  loss  through  the  air 
boundary  layer.  Transport  though  the  air  boundary 
layer  is  controlled  by  diffusion  and  limited  by  the 
transfer  of  chemical  from  the  aqueous  phase  to  the 
gas  phase  by  the  very  low  Kh- 


Depth  (cm) 

Figure  8.  Upward  Transport  and  Development  of 
a  Surface  Layer 

Figure  10  shows  the  effect  of  developing  the 
enhanced  surface  layer  with  60  days  of  evaporation 
(>0.5  cm/day),  followed  by  precipitation  for  5  days 
(0.5  cm/day).  The  enhanced  surface  layer  found  in 
the  top  0.1  cm  of  soil  is  transported  down  leaving  just 
a  small  enhancement  at  a  depth  of  about  0,5  cm. 
Another  simulation  was  run  that  included  the  same 
evaporation  and  precipitation,  but  was  followed  by 
another  5  day  evaporation  period  (-0.5  cm/day)  and 
the  surface  enhancement  returned  at  about  the  same 
concentration. 


Figure  9.  Detail  of  the  Surface  Layer  Formed  by 
Evaporation 


SoM  Concentration  (ug/cm3) 


Figure  10.  Effect  of  Precipitation  on  the 
Enhanced  Surface  Layer 


DISCUSSION 

As  a  result  of  these  initial  results,  a  more  detailed 
mechanistic  numerical  model  is  being  developed. 
This  model  is  being  based  on  TOUGH2  (Pruess, 
1991)  with  modifications  pertinent  to  the  landmine 
application  and  will  be  called  T2TNT.  Modifications 
being  made  or  planned  to  be  made  include: 

1 .  Addition  of  TNT  and  DNT  vapor  components.  - 
Landmines  typically  emit  TNT  and  DNT  vapors. 
The  vapor  pressure  of  DNT  is  higher  than  TNT 
and,  if  present,  will  probably  reach  any  chemical 
sensor  before  TNT. 

2.  Dusty  Gas  Model  for  gas  diffusion.  -  Gas 
diffusion  is  a  dominant  transport  mode  for  TNT 
and  DNT  vapors  in  the  subsurface.  Therefore, 
the  Dusty  Gas  Model  (Webb,  1998)  will  be 
implemented. 

3.  Liquid  diffusion  of  dissolved  TNT  and  DNT.  - 
Liquid  diffusion  is  not  present  in  the  standard 
version,  although  some  special  EOS  modules 
include  it.  Liquid  diffusion  analogous  to  gas 
diffusion  will  be  included  because  of  the 
significant  chemical  concentration  in  the  liquid 
phase. 

4.  Partition  coefficient  as  a  function  of  saturation.  -- 
The  solid  partition  coefficient  may  be  a  strong 
function  of  saturation,  especially  at  low  moisture 
content  where  the  partition  coefficient  may 
increase  dramatically  (Petersen,  et  al,  1995). 

5.  Boundary  layer  specifications  for  transport  at 
soil  surface.  -  For  transport  from  the  soil  to  the 
atmosphere,  a  boundary  layer  for  heat  and  mass 
transfer  will  be  implemented.  The  exact 
modifications  are  still  being  developed. 

6.  Precipitation  and  evaporation  boundary 
conditions-  -  Precipitation  boundary  conditions 
will  be  added.  The  evaporation  boundary 
condition  may  simply  involve  specification  of  a 


boundary  layer  and  a  boundary  relative  humidity, 
which  will  be  time  dependent. 

7.  Diurnal  and  seasonal  variations  in  atmospheric 
conditions  -  In  order  to  simulate  daily  and 
seasonal  fluctuations,  time-dependent  boundary 
conditions  for  the  pressure,  temperature,  and 
relative  humidity  will  be  implemented. 

The  resulting  code  will  be  used  to  develop  an 
effective  operational  strategy  for  the  design  and 
deployment  of  landmine  chemical  sensors.  The  code 
will  also  be  used  within  the  ITOUGH  framework  to 
assist  in  the  design  of  column  experiments  to  be 
conducted  at  New  Mexico  Tech  during  the  next  few 
years.  T2TNT  will  play  an  important  part  in  the 
effective  use  of  chemical  sensors  for  landmine 
detection  and  removal. 

CONCLUSIONS 

The  environmental  fate  and  transport  of  chemical 
signatures  from  landmines  is  important  for  the  design 
and  operation  of  chemical  sensors.  The  explosive 
vapors  are  predominantly  found  sorbed  to  soil 
particles  or  in  the  liquid  phase;  only  a  small  fraction 
is  present  in  the  gas  phase.  As  a  result,  diffusion  and 
advection  of  the  liquid  water  dominates  the  transport 
of  the  chemicals  from  the  buried  landmine  to  the  land 
surface.  Precipitation  and  evaporation  also  strongly 
influence  the  movement  of  the  chemical  signature. 

The  results  of  initial  screening  studies  have 
confirmed  the  influence  of  environmental  conditions 
and  soil  parameters.  The  burial  depth  of  the 
landmine  is  a  significant  factor.  For  shallow 
landmines,  the  chemical  appears  much  sooner  and  at 
a  much  higher  concentration  than  for  deeper 
landmines.  Precipitation  and  evaporation  have  a 
significant  effect  on  the  transport  of  TNT  in  the 
subsurface.  The  chemical  concentration  at  the 
surface  varies  by  many  orders  of  magnitudes 
depending  upon  whether  precipitation  or  evaporation 
is  occurring.  Under  evaporation  conditions,  a  surface 
“crust”  can  form  where  the  surface  concentration  is 
higher  than  in  the  subsurface.  These  conditions 
would  greatly  enhance  the  detection  capability  of 
chemical  sensors.  As  a  result,  a  mechanistic 
numerical  model  based  on  TOUGH2  (Pruess,  1991), 
called  T2TNT,  is  currently  being  developed.  T2TNT 
will  include  a  number  of  modifications  and 
enhancements  that  should  be  of  general  interest  to 
many  TOUGH2  users.  The  use  of  T2TNT  in 
understanding  the  environmental  fate  and  transport  of 
TNT  in  the  soil  will  contribute  to  the  improved 
design,  performance  and  operation  of  chemical 
sensors  in  the  detection  of  buried  landmines. 


NOMENCLATURE 
a  air  volume  fraction  =  0  -  0 
D  diffusion  coefficient 
IQ  distribution  coefficient  for  sorption 
Kh  Henry’s  Law  constant 
0  volumetric  moisture  content 
(()  porosity 
p  density 

Subscripts 
E  effective 
g  gas 
1  liquid 

Superscripts 
a  air 
w  water 
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One  Dimensional  Transport  of  Vapor  From 
a  Buried  Land  Mine 

William  Jury  and  Lei  Guo,  University  of  California,  Riverside 

Introduction 

Constructing  a  general  strategy  for  detecting  emissions  at  the  soil  surface  from  buried  land 
mines  requires  an  understanding  of  the  effect  of  soil  and  environmental  conditions  on  signal 
strength  and  location.  While  too  simple  to  simulate  in  detail  the  actual  transport  aind  reaction 
processes  occurring  in  a  natural  setting,  screening  models  are  useful  tools  for  elucidating  the 
relationships  between  soil  parameters  and  output  characteristics  of  interest.  This  first  model 
is  one  dimensional,  representing  steady  emissions  from  a  buried  source  of  TNT  vapor.  It  is  an 
extension  of  the  screening  model  papers  of  Jury  et  al.  (1983;  1984abc;  1990) 

Model  Description 


Transport  Equations:  The  analysis  of  the  original  screening  model  paper  (Jury  et  al.  1983) 
produces  the  familiar  transport  equation 


dCt 

dt 


=  Db 


d^Ct 

dz^ 


Ve 


dCt 

dz 


—  ^eCi 


(1) 


where 

Ci  =  pbCa  +  OCc  +  CiCg  (2) 

is  the  total  concentration  of  TNT  in  the  soil,  Ca^Cc^Cg  are  the  concentrations  in  the  sorbed, 
dissolved,  and  vapor  phases,  respectively,  pb  is  soil  bulk  density,  and  6  is  volumetric  water 
content. 

The  effective  diffusion  coefficient  De  is  the  combined  diffusion  in  the  dissolved  and  gaseous 
phases,  given  by 

(3) 

where  Kf{  is  the  dimensionless  Hemy’s  constzmt,  Dg''',  Df  axe  the  binary  diffusion  coefficients 
of  the  chemical  in  air  and  water  respectively,  and  ^5,  ^,1,  are  the  tortuosity  factors,  given  by  the 
modified  Millington-Quirk  tortuosity  law 


(Jin  and  Jury,  1993).  The  effective  solute  velocity  is  given  by 


Vf  =  — 
^  R 


where  J^)  is  the  water  flux  rate. 

The  retardation  factor  R  is  given  by 


i?  =  1  + 


pbKd 

e 


(4) 

(5) 


(6) 


where  Kd  is  the  slope  of  the  linearized  adsorption  isotherm  (the  distribution  coefficient). 


The  upper  boundai-y  condition  is  assumed  to  be  controlled  by  diffusion  through  a  stagnant 
air  boundary  layer  above  the  soil  surface,  and  is  expressed  as: 


-DB^  +  VECt  =  HECt 
az 


2  =  0 


where  He  is  the  effective  mass  transfer  coefficient  through  the  boundary  layer,  given  by 


(7) 


He  = 


KHDf 

JRdbi 


(8) 


where  is  the  boundary  layer  thickness,  typically  about  0.5  cm  (Jury  et  al.  19S4c). 

Constant  Flux  Emission  From  the  Land  Mine:  The  lower  boundary  condition  describing 
emission  of  a  constant  flux  from  the  land  mine  is  expressed  mathematically  as 


-De^  +  VeCi  =  F  z=^-P  (9) 

oz 

where  F  is  the  flux  of  gas  across  the  plane  at  2  =  — P. 

The  system  described  by  these  equations  will  come  to  a  steady  state  eventually.  We  can 
easily  solve  for  the  steady  profile  and  flux  emission  from  the  surface.  At  steady  state,  the 
concentration  at  2  =  ~P  will  also  be  constant,  so  that  the  solution  will  also  apply  to  the  case 
where  the  bottom  of  the  profile  is  held  at  a  constant  value. 

The  steady  form  of  the  transport  equation  is 


which  has  the  solution 


C,{z)  =  A  e.xp[^(l  -01  +  5  exp[ 


'2De' 


2De 


d  +  O]  ;  + 


(11) 


The  corresponding  general  solution  for  the  flux  is  given  by 

Js(^)  =  Y  [-4(1  +  0  e>^p[^(l  -  01  +  -5(1  -  0  +  Ol]  (12) 

The  boundary  conditions  (2-3)  can  be  applied  to  (11)-(12),  producing 
Ct(2)  =  ^  [{2He  +  V£(e  -  1))  exp[n2(l  -  0]  -  {2He  -  Ve{^  +  1))  exp[n2(l  +  0]]  ^  (13) 

M^)  =  I  [(277^(1  +  ?)  +  ^^^)  exp[n2(l  -  0]  +  i2HE{^  -  1)  -  ^^^)  exp[n2(l  +  0)]  (14) 

where 

13  =  i2HE{l  +  0  +  ^^^)  exp[-nP(l  -  01  +  (2f7e(e  -  1)  -  ^^^)  exp[-nP(l  +  0]  (15) 

Of  particular  interest  are  the  flux  from  the  surface 


Js{0)  = 


AHe^F 

P  . 


(16) 


and  the  concentration  at  z  =  —  P 


Ct{-P)  =  Co  = 


2FTj 

PVe 


(17) 


2 


where 


Tj  =  (2Hb  +  Vb{^-1))  exp[-fiP(l  -  0]  -  {2Hb  -  V£:(^  +  1))  exp[-J7P(l  +  0]  (IS) 

Note  that  the  surface  concentration  is  given  by  Js(,0)IHb- 

We  can  rewrite  (16)  in  terms  of  three  dimensionless  variables 

_ I _  (19) 

~  F  [l+l/^+2J/(w7)]e-T'/2*(i-«)  +  [l-l/^-2(J/(a;7)]e-''/2‘(i+«) 


where 


PVb 


In  the  case  where  there  is  no  evaporation  (Vg  =  0),  this  reduces  to 


(20) 


_  Js(0)  _ _ 2 _  k  =  p/K  e  = 

F  (1  +  Ke)e'‘  +  (1  —  /ce)e  ''  V  Dg  PHe 


(21) 


Focusing  first  on  (21),  we  interpret  k  as  the  ratio  of  the  diffusion  time  (i.e.  the  travel  time 
to  the  surface)  to  the  decay  time.  A  large  value  of  k  will  mean  that  most  of  the  flux  will  be 
consruned  en  route  to  the  surface.  The  second  dimensioirless  variable  c  can  be  interpreted  as 
the  ratio  of  the  transfer  coefficient  to  the  boundary  layer  and  the  transfer  coefficient  through 
the  boundary  layer.  Using  all  of  the  definitions  for  these  coefficients,  we  can  express  e  as 


<  =  +  «.(«))  <  = 


(22) 


where  ^v(o)  Is  the  tortuosity  factor  for  the  gas  phcise  in  soil,  ^w(^)  is  the  tortuosity  factor 
for  the  dissolved  phaise  in  soil,  dg  ~  0.5  cm  is  the  thickness  of  the  air  boundary  layer J^//  is 
Henrys  constant,  and  D"”"  are  the  diffusion  coefficient  of  the  chemical  in  water  and  air, 

respectively.  The  parameter  C  expresses  the  chemical’s  preference  for  the  dissolved  or  vapor 
phase. 

Constant  Concentration  at  the  Depth  of  the  Land  Mine:  The  boundary  condition  for 
the  case  of  a  constant  concentration  at  the  lower  boundaiy  is  given  by 


Ct(-P)  =  Co  (23) 

We  may  write  down  the  solutions  for  the  steady  profile  and  flux  for  the  case  of  constant 
concentration  (23)  using  (17) 

Ct(z)  =  ^  [{2He  +  Vb{^  -  1))  exp[n2(l  -  01  -  (2Pg  -  Vb{^  +  1))  exp[flz(l  +  Oj]  (24) 
Js{z)  =  ^  [(2P£(1  +  0  +  ^^^)  exp(f22(l  -  01  +  {2He{^  -  1)  -  ^^^^)  exp[flz(l  +  Ol]  (25) 

7^(0)  =  (26) 

This  solution  will  later  be  superimposed  on  the  constant  flux  solution  to  allow  an  analysis 
of  the  relative  importance  of  the  two  signals  for  typical  values  of  the  boundary  parameters. 
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Results 


Constant  Flux  Emissions:  Calculations  are  performed  for  TNT  using  values  of  Kh  = 
6  X  10“^  and  Koc  =  300  cm^/g  for  a  range  of  half  lives.  Figure  1  shows  the  predicted  relative 
flux  from  the  surface  as  a  function  of  the  steady  upward  evaporation  rate,  showing  a  significant 
enhancement  due  to  upwai'd  advection.  Note  that  the  rate  of  loss  is  higher  when  the  water 


Evaporation  Rate  (cnVd) 


Figure  1:  Flux  from  the  surface  as  a  function  of  evaporation  rate  for  the  case  of  a  buried  mine  at 
10  cm  depth.  Soil  conditions:  pi  =1.5  g/cm^,  foe  =  -01  Chemical  properties;  itroc=300  cm^/g, 
=  10-^  ti/2  =  1  yr-i 

content  is  lower.  This  is  because  the  transport  is  primarily  in  the  dissolved  phase,  and  the 

travel  time  to  the  surface  by  advection  is  inversely  proportional  to  6. 

Figure  2  shows  the  predicted  relative  flux  from  the  surface  as  a  function  of  water  content 
for  various  values  of  evaporation.  At  high  evaporation  rate,  the  loss  is  dominated  by  advec- 


Figure  2:  Flux  from  the  surface  as  a  function  of  water  content  for  the  case  of  a  buried  mine  at 
10  cm  depth.  Soil  conditions;  pb  =1.5  g/cm^,  foe  =  -01  Chemical  properties:  i<'oc=300  cm^/g, 
10-^  ri/2  =  1  yi-i 

tion,  and  therefore  increases  with  decreasing  6.  At  low  evaporation,  the  system  shifts  over  to 
diffusion-dominated,  showing  a  modest  increase  with  increasing  6.  Vapor  diffusion  is  insignifi¬ 
cant  compared  to  liquid  diffusion  over  much  of  the  range  of  water  content  because  of  the  low 
Henry’s  constant.  The  ratio  r  of  liquid  to  vapor  diffusive  transport  for  a  given  concentration 
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gradient  is  given  by 


Dsoit  _  ^  10-“,  5  2  c-i-i  ^ 

''~K}iDl°‘‘  KuDf%(a)~  Kii\l-Sy  4>  <l> 

Table  1  shows  various  values  of  r  calculated  with  (27)  as  a  function  of  S. 


s 

r 

S 

r 

1.0 

oo 

0.2 

10 

0.6 

375 

0.071 

1 

0.4 

74 

0.0 

0.0 

Table  1:  Ratio  of  liquid  to  vapor  diffusion  for  TNT  as  a  function  of  water  saturation 

Figure  3  shows  the  predicted  relative  flux  from  the  surface  as  a  function  of  organic  carbon 
fraction. 


Figui'e  3:  Flux  from  the  surface  as  a  function  of  organic  carbon  fraction  for  the  case  of  a  buried 
mine  at  10  cm  depth.  Soil  conditions:  =1.5  g/cm^,  6  =  0.3;  Chemical  properties:  ifoc^^OO 

cm^/g,  =  6  X  10-^  Ti/2  =  1  yr"^ 

Since  the  value  used  for  Koc  Is  300,  the  range  of  Kd  simulated  in  this  figure  is  0  to  15.  This 
figure  holds  no  surprises,  merely  illustrating  the  effect  of  increasing  residence  time  on  decreasing 
volatilization  loss  because  of  degx'adation. 

Figure  4  shows  the  predicted  relative  flux  from  the  surface  as  a  function  of  the  chemical 
half  life.  Half  life  is  obviously  a  sensitive  parameter,  and  unfortunately  is  also  one  for  which 
no  data  is  available.  Significant  in  Figure  4  is  the  much  smaller  loss  when  no  upward  flow 
is  occurring.  The  enhanced  loss  under  upward  water  flow  is  caused  by  accumulation  of  TNT 
at  the  surface  by  deposition  from  the  evaporation  process,  which  causes  a  buildup  of  vapor 
concentration  and  an  increctsed  driving  force  for  vapor  diffusion  across  the  boundary  layer. 
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Figure  4:  Flux  from  the  surface  as  a  function  of  chemical  half  life. 

Figure  5  shows  the  predicted  relative  flux  from  the  surface  as  a  function  of  water  content 
for  the  case  of  no  evaporation  (note  expanded  scale). 


Figure  5:  Flux  from  the  surface  as  a  function  of  water  content  and  no  evaporation. 

Figure  6  shows  the  falloff  in  flux  with  depth  of  land  mine,  which  is  not  too  significant  in 
the  5-15  cm  range. 


Figure  6:  Flux  from  surface  as  a  function  of  land  mine  depth. 

The  Surface  Concentration:  According  to  the  model,  the  total  surface  concentration  is 
given  by  Js{0)IHe-  Thus,  all  of  the  graph  results  for  the  surface  flux  are  proportional  to  the 
surface  conceiitration.  Since  the  Y-axis  gives  relative  flux,  we  can  convert  to  actual  surface 
concentration  by  multiplying  the  axis  value  by  F/He  As  a  rough  estimate,  we  can  use  the 
following  data  (Table  2)  from  the  report. 
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Parameter 

Value 

Units 

Kh 

6  X  10-  ' 

> 

Ka 

5 

cm^/ g 

Re 

15 

- 

h 

8600 

cm/d 

He 

3.5  X  lO"'' 

cm/d 

F 

10-16  _  10-18 

g/cm^/s 

FIHe 

0.25-25 

Mg/L 

Table  2:  Representative  parameters  for  TNT  adopted  from  the  report:  SAND97-14^6 


The  concentration  profile  in  the  soil  can  build  up  substantially  near  the  surface  because  the 
boundary  layer  is  the  limiting  step  to  loss  from  the  soil.  Figure  7  shows  a  plot  of  concentration 
(log  scale)  versus  depth  as  a  function  of  water  content  for  the  highest  rate  (E=0,5  cm/d)  of 
evaporation. 


Figure  7:  Soil  concentration  profile  as  a  function  of  water  content  for  the  case  of  a  buried  mine 
at  10  cm  depth.  Soil  conditions:  E  =  0.5  cm/d,  pb  =1.5  g/cm^,  foe  =  0.01;  Chemical  properties: 
i^oc=300  cm^/g,  Kh  =  6  X  10”^ 

As  the  evaporation  rate  declines,  the  profile  becomes  less  extreme,  as  shown  in  Figure  8. 

Discussion 

The  one  dimensional  simulations  are  of  limited  value  for  estimating  the  magnitude  of  the  surface 
flux,  but  they  do  show  the  main  features  the  govern  the  true  process.  Compounds  such  as  TNT, 
with  extremely  low  Heniys  constant,  will  be  substantially  influenced  by  upward  flow,  both 
because  advection  will  provide  the  most  rapid  transport  pathway  to  the  surface,  and  because 
the  air  boundary  layer  will  impede  removal  from  the  surface. 

The  loss  of  chemical  under  zero  evaporation  forms  a  lower  limit  to  accumulation  at  the 
surface,  where  over  time  the  concentration  in  the  soil  will  form  a  decreasing  gradient  to  the 
surface.  Table  3  shows  various  values  of  HeC{0)/F  for  the  case  of  zero  evaporation,  a  water 
content  of  0.3,  a  depth  of  10  cm  for  the  upper  surface  of  the  land  mine,  and  an  organic  carbon 
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Figure  8:  Soil  concentration  profile  as  a  function  of  evaporation  rate  for  the  case  of  a  biuried 
mine  at  10  cm  depth.  Soil  conditions:  6  =0.3,  pb  =1.5  g/cm^,  foe  =  0.01;  Chemical  properties: 
i(:oc=300  cm^/g,  Kh  =  Qx  10-^ 

fraction  of  0.01.  These  can  be  multiplied  by  the  values  of  F/He  in  Table  to  estimate  surface 
concentrations. 


T 

HeC{0)/F 

3 

.072 

2 

.045 

1 

.014 

.5 

.003 

.25 

.0003 

Table  3:  Estimated  values  of  the  dimensionless  surface  concentration  for  various  values  of  the 
half  life  for  the  case  of  zero  upward  flow  of  water  to  the  surface. 


Constant  Concentration  at  Lower  Boundary:  The  relation  between  the  flux  value  and 
the  corresponding  concentration,  as  expressed  in  (17),  is  given  in  Figure  9  for  representative 
conditions. 

This  relationship  is  nearly  linear  under  a  range  of  soil  and  environmental  conditions.  For 
the  case  where  evaporation  is  nonzero,  it  becomes 


[HeVe{1  +  ^)  +  ‘^.peDe]  ^ 
~  l2HE  +  VEi^-l)] 


(28) 


and  if  evaporation  is  zero  it  becomes  simply 

F  w  '/peF^Cq  (29) 

For  the  conditions  represented  in  Figui'e  9,  the  range  of  F  values  (10“^®  —  10“^®  g/cm^/s) 
typical  of  land  mine  emissions  correspond  to  soil  concentrations  Co  at  the  lower  boundary  of 
7.1  X  10“'^  —  7.5  X  10“®  p  g/cm^.  In  other  words,  a  soil  concentration  of  this  magnitude  would 
give  rise  to  the  same  surface  emissions  as  the  flux  in  this  range.  Figure  10  illustrates  the  case 
where  both  land  mine  flux  emission  of  TNT  and  constant  soil  concentrations  are  present,  where 
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Figiire  9:  Relation  between  the  flux  F  from  depth  z  —  —P  and  the  concentration  Co  at  that 
depth.  Solid  vei'tical  line  brackets  the  range  of  fluxes  F  «  10“^®  — g/cm^/s  given  in  report 
SAND97-1426  .  Soil  conditions:  pi  =1.5  g/cm^,  6  =  0.3,  foe  =  0.01;  Chemical  properties: 
i!roc=300  cm^/g,  Kh  =  6  X  10“^  ti/2  =  1  yr~^ 


F=  10“^*^  g/cm^/s  and  Co  =  4.6  x  10“^  p  g/cm^,  the  latter  taken  from  the  screening  calculations 
of  the  Samdia  report  SAND97-1426  .  As  seen  in  this  flgure,  under  these  conditions  the  signal 
from  the  soil  concentrations  at  10  cm  would  dominate  the  surface  emissions. 


Figm'e  10:  Flux  at  the  surface  arising  from  both  land  mine  flux  emissions  F=  10  g/cm^/s 
and  soil  concentrations  Co  =  4.6  x  10~^  p  g/cm^  at  z  =  —P. 
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Transient  Conditions:  The  transient  problem  is  more  complex  to  analyse,  but  an  approxi¬ 
mate  estimate  of  the  time  required  to  reach  equilibrium  can  be  calculated  for  the  case  of  zero 
evaporation  (the  longest  time  to  reach  equilibrium)  from  the  diffusion  time 


16I>£ 


Table  4  gives  the  equilibration  times  calculated  for  TNT  using  Eq.  18. 


(30) 


e 

td  (d) 

.1 

919 

.2 

244 

.3 

115 

.4 

68 

.5 

46 

Table  4:  Approximate  time  to  reach  steady  state  for  the  case  of  zero  evaporation  for  the 
conditions  foe  =  0.005,  P  =  10.,  K}j  =  6  x  10“^,  Koc  =  300. 


Prior  to  these  times  after  insertion  of  the  signal  at  the  depth  of  the  land  mine,  the  fluxes 
and  concentrations  at  the  surface  would  be  less  than  the  amounts  calculated  from  the  steady 
model. 
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Two  Dimensional  Transport  of  Vapor  from 
a  Buried  Land  Mine 

William  Jury  and  Lei  Guo,  University  of  California,  Riverside 

Introduction 

Constructing  a  general  strategy  for  detecting  emissions  at  the  soil  surface  from  buried  land  mines 
requires  an  understanding  of  the  effect  of  soil  and  environmental  conditions  on  signal  strength 
and  location.  While  too  simple  to  simulate  in  detail  the  actual  transport  and  reaction  processes 
occurring  in  a  natural  setting,  screening  models  are  useful  tools  for  elucidating  the  relationships 
between  soil  pai'ameters  and  output  characteinstics  of  interest.  In  a  previous  report  we  examined 
the  release  of  TNT  vapor  from  a  subsurface  land  mine  using  a  one-dimensional  model.  In  this 
report  we  use  a  two  dimensional  model,  representing  steady  emissions  from  a  buried  plane  source 
of  TNT  vapor. 

Figure  1  shows  a  schematic  of  the  system  represented  with  the  model.  The  land  mine  is  character¬ 
ized  by  a  constant  flux  plane  at  ^  =  -P,  which  emits  a  constant  upward  flux  within  —W  <  x  <W, 
and  zero  upward  flux  outside  of  this  range.  At  a  substantial  lateral  distance  x  =  ±L,  the  hori¬ 
zontal  flux  is  set  equal  to  zero.  At  the  surface,  chemical  vapor  can  volatilize  by  diffusing  through 
a  stagnant  air  boundary  layer  of  thickness  Sz  =  above  which  the  air  is  assumed  to  be  at  zero 
chemical  concentration. 


2L 


2  W 


Figure  1:  Schematic  of  the  soil  regime  represented  by  the  model 


Transport  Equations 


Two  dimensional  transport  of  TNT  by  vapor  and  liquid  diffusion  and  liquid  advection  may  be 
described  by 


dt 


^  d-c,  ^  ^Ct  ,,  dc, 


(1) 


where 

Ct  =  PbCa  +  OCc  +  dCg 


(2) 


is  the  total  concentration  of  TNT  in  the  soil,  Ca,Cc,Cg  are  the  concentrations  in  the  sorbed, 
dissolved,  and  vapor  phases,  respectively,  pt  is  soil  bulk  density,  and  6  is  volumetric  water  content. 

The  effective  diffusion  coefficient  Pe  is  the  combined  diffusion  in  the  dissolved  and  gaseous 
phases,  given  by 


De  = 


R 


(3) 


where  Kh  is  the  dimensionless  Henry’s  constant,  Dg",  Df  are  the  binary  diffusion  coefficients  of  the 
chemical  in  air  and  water  respectively,  and  ^gi^w  th®  tortuosity  factors,  given  by  the  modified 
Millington-Quirk  tortuosity  law 


^(1)  = 


(/.4/3’ 


X  =  6,a 


(4) 


(Jin  and  Jury,  1993). 

Figure  2  shows  a  plot  of  Pe  versus  9  for  TNT  for  several  values  of  organic  carbon  fraction  foe- 
Because  of  its  low  Henrys  constant,  TNT  is  dominated  by  liquid  diffusion  over  most  of  the  range 
of  water  content,  and  has  an  effective  value  that  is  nearly  200  times  larger  at  saturation  than  at 
extremely  low  water  content. 


TJ 


Figure  2:  Effective  diffusion  coefficient  of  TNT  calculated  using  Eq.  (3). 

The  effective  solute  velocity  is  given  by 

Ve  =  — 

^  R 

where  is  the  water  flux  rate,  assumed  to  be  one  dimensional  in  the  vertical  direction. 
The  retardation  factor  R.  is  given  by 


(5) 


R  —  1  + 


PbKd 


(6) 


2 


where  Kd  is  the  slope  of  the  Hirearized  adsorption  isotherm  (the  distribution  coefficient). 

The  upper  boundaiy  condition  is  assumed  to  be  controlled  by  diffusion  through  a  stagnant  air 
boundary  layer  above  the  soil  surface,  and  is  expressed  as. 


-DE^+VECt  =  HECt 
OZ 


z  =  0;  0  <  rr  <  L 


where  He  is  the  effective  mass  transfer  coefficient  through  the  boundary  layer,  given  by 


He  = 


khd: 


atr 
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Rdbi 


(7) 


(8) 


where  dy  is  the  boundary  layer  thickness,  typically  about  0.5  cm  (Jury  et  al.  1984c). 

Constant  Flux  Emission  From  the  Land  Mine:  The  lower  boundary  condition  describing 
emission  of  a  constant  flux  from  the  land  mine  is  expressed  mathematically  as 


-DE^  +  VECt  =  f{x)  z  =  -P;Q<x<L  (9) 

OZ 

where  f(x)  is  an  arbitraiy  function  representing  the  upward  flux  of  TNT  across  the  plane  at  z  —  —P. 

The  planes  formed  by  the  lines  x  =  0  and  i  =  L  are  planes  of  symmetry,  which  obey  the 
conditions 

dCt 
dx 

dCt 

dx 

The  system  described  by  these  equations  will  come  to  a  steady  state  eventually.  We  can  solve 
for  the  steady  profile  and  flux  emission  by  using  separation  of  variables  and  Fourier  series.  We  first 

let 


=  0 

a;  =  0; 

-P  <  z  <0 

(10) 

=  0 

x  =  L; 

-P  <  z  <  0 

(11) 

Ctix.z)  =  Z{z)Xix) 

and  plug  this  into  the  steady  state  form  of  the  transport  equation,  with  the  result 

Rearranging,  we  obtain 

DEd?X  VEdZ  DEfZ__  2 

X  dx^-  Z  dz  Z  dz'^ 

where  a;  must  be  a  constant.  This  produces  two  equations  to  solve. 

Db^  +  ij‘^x  =  a 

Db^-Ve^-(<^^  +  I‘e)Z  =  0 
The  general  solution  to  (15)  is 

X{x)  =  Ax  cos(/0x)  +  Bx  cos(/3x)  ;  d  = 


(12) 

(13) 

(14) 

(15) 

(16) 

(17) 
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and  the  general  solution  to  (16)  is 


The  boundary  conditions  (5-6)  can  be  applied  to  (17),  producing 


(IS) 


Ax=0  and  P  — =  0,1,2, 


(19) 


This  result  mezms  that  there  are  an  infinite  number  of  values  of  ^  which  are  functions  now  of 
N.  Subsequently,  any  parameter  depending  on  N  will  have  a  subscript.  In  particular,  all  Xa'(x) 
satisfying 

Xn{x)  =  BxN  cosi^)  N  =  0,1,2,...  (20) 

obey  both  the  equation  and  the  boundary  conditions  governing  X(x). 

Applying  the  condition  (7)  to  (18),  we  obtain 


[^(1  +  M  -  Hb]A.  +  [^(1  -  Cn)  -  He]B,  =  0  -4  =  jnAx 


where 


7N  = 


2He-Ve{1^^n) 


(21) 


(22) 


I4;(l  -  ^n)  -  277£ 

and  the  most  general  form  of  Z{z)  satisfying  the  equation  and  the  surface  boundary  condition 


is 


Z^r{z)  =  A;A»{exp[^^(l  -  ^aO]  +  IN 


(23) 


By  linear  superposition,  the  sum  of  solutions  to  the  linear  transport  equation  is  also  a  solution. 
Therefore,  the  most  general  solution  to  the  original  equation  and  all  of  the  boundary  conditions 
used  thus  far  is 


Ct  (x,z)  =  E^^o■^^rZ^r  (z)XAr(x)  Ef,i  =  A^nBxN 


(24) 


where 


Jwz{x,z)  =  T,'^,-QE^TQf.{{z)X^J{x) 

Qn{z)  =  ^[(1  +  (n)  w(l  -  ^n) 

Plugging  in  for  the  functions,  we  have 

Ct(x,z)  =  ~  ^ 


(25) 

(26) 

(27) 


Equation  (27)  has  an  infinite  number  of  unknowns  which  are  evaluated  by  applying  the 
lower  boundary  condition  (9)  to  the  equation  and  making  use  of  the  orthogonal  property  of  the 
cosine  functions.  Thus, 


.Nttx, 

fix)  =  0-Ea'Qa'(--P)cos[-^] 


(28) 
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where  Qr,r{-P)  is  calculated  with  (20).  Equation  (28)  is  a  Fourier  cosine  series,  and  the  coefficients 
En  can  be  evaluated  by  multiplying  both  sides  by  cos(M7rx/L)  and  integrating  from  0  to  L.  This 
yields 


Eo 


1 

LQoi-P) 


(29) 


(30) 

LQn{-P)Jo  -  L  ’ 

It  is  also  useful  to  calculate  the  flux  at  the  surface  as  a  function  of  position.  To  do  this  we  plug 
(27)  into  the  expression  for  the  flux  and  evaluate  it  at  z=0.  This  produces 


J(x,0))  =  -De^U=o  +  ^£^(0)  =  i:^^oENXN(x)n]^  Qn  =  ^[(l  +  (n)  +  iNil  -  ^n)]  (31) 

Constant  Concentration  at  the  Depth  of  the  Land  Mine:  Unlike  the  one  dimensional 
case,  the  constant  flux  and  constant  concentration  solutions  are  not  mere  multiples  of  each  other. 
In  the  constant  flux  case  above,  for  example,  the  concentration  along  0<a:<W^atz  =  — Pis  not 
constant.  Thus,  the  constant  concentration  problem  must  be  solved  separately.  The  lower  bound¬ 
ary  condition  describing  a  constant  concentration  along  the  plane  of  the  land  mine  is  expressed 
mathematically  as 

Ct=-9{x)  z^--P]0<x<L  (32) 

where  9{x)  is  an  arbitrary  function  representing  the  chemical  concentration  across  the  plane  at 
2  =  -P. 

The  rest  of  the  boundary  conditions  are  as  before,  and  the  method  of  solution  is  again  separation 
of  variables,  producing  the  same  genei'al  solution  (27).  The  unknowns  Eisj  are  evaluated  in  this  case 
by  applying  the  lower  boundary  condition  (32)  to  the  equation  and  making  use  of  the  orthogonal 
property  of  the  cosine  functions.  Thus, 


(7(x)-S~o^Ar^A'(-^)cos[  ^  ] 

(33) 

where  t/  T/  P 

Zwi-P)  -  exp[^^  l)]  +  7Nexp[  (1  +  Cn)]} 

(34) 

Equation  (33)  is  a  Fourier  cosine  series,  and  the  coefficients  can  be  evaluated  by  multiplying 

both  sides  by  cos{MTTxf L)  and  integi'ating  from  0  to  1/  as  before.  This  yields 

=  LZ«(-P)  i 

(35) 

The  flux  at  the  surface  is  simply  equal  to  iy£;Cf(x,0),  or 

(36) 

PI  TTX 

J(x,0))-i7£E~o^Ar(l  +  7N)cos[  ^  j 

(37) 

We  may  now  solve  the  problem  for  particular  choices  for  the  flux  f{x)  or 
across  the  boundary  at  2  =  ~P. 

concentration  p(x) 
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Results 


Plane  Source  Flux  Emission:  In  this  case 


fix)  =  F 
fix)  =  0 


Applying  this  to  (29)-(30),  we  obtain 


Eq  = 


0  <  X  <  VF 

W  <x<L 


FW 


LQo 


= 


2F 


A'ttQat 


N-nW 

sm(^^)  AT  =  1,2,.. 


and  thus  that 


,  FWZq{z)  .  2F^^  ,N'!tx^  .  ^N'KW^ZMiz) 


J(x,0))  = 


FWU,  .  2F^^  cosF^’^^’ 

■“N-X  COS[ 


ATTrW,  VLn 


H - cos[  -  ]sin[  ] 


(38) 

(39) 

(40) 

(41) 

(42) 


TQo  '  L  L  ^NQ^^ 

Figure  3  shows  a  plot  of  the  relative  flux  J„,r(x,0)/F  as  a  function  of  position  on  the  surface  for 
the  extreme  case  of  infinite  half  life  for  TNT,  high  water  content  (creating  the  maximum  diffusion 
coefficient),  and  high  evaporation  rate. 


Figure  3:  Flux  from  the  surface  as  a  function  of  lateral  position  for  the  Ccise  of  a  buried  mine  of 
width  2W  =  30  cm  at  10  cm  depth. 

Figure  4  shows  the  effect  of  decreasing  evaporation  rate  on  the  lateral  spreading.  As  E  declines, 
more  lateral  movement  takes  place  because  the  travel  time  to  the  surface  increases.  This  figure  is 
possibly  misleading  however,  because  of  the  assumption  of  zero  degradation.  Whenever  degradation 
is  present,  chemicals  following  paths  with  long  travel  times  (as  are  the  paths  lateral  to  the  vertical 
plane  formed  by  the  boundaries  of  the  land  mine)  are  attenuated  prior  to  arrival 

Figure  5  is  a  repeat  of  Figure  4,  except  that  the  half  life  of  TNT  is  assumed  to  be  1  year. 
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Figure  4;  Effect  of  evaporation  rate  on  the  flux  from  the  surface  as  a  function  of  lateral  position. 
Same  conditions  as  in  Fig.  3 


Figm-e  5:  Effect  of  evaporation  rate  on  the  flux  from  the  surface  as  a  function  of  lateral  position. 
Same  conditions  as  in  Fig.  4  except  that  r  =  1  yr. 

In  contrast  to  the  case  where  no  degradation  is  occurring,  the  lateral  extent  of  the  surface  flux 
is  much  less. 

Water  content  also  has  a  significant  effect  on  signal  strength  and  lateral  spreading.  Figure  6  is 
a  recalculation  of  Figure  5  with  the  water  content  reduced  from  0.5  to  0.25. 

In  this  case,  the  advective  travel  time  to  the  surface  (which  is  inversely  proportional  to  6)  is 
decreased,  causing  higher  surface  flux  above  the  land  mine  when  E  ^  0.  At  the  same  time,  however, 
diffusive  travel  time  is  lengthened  significantly  (see  Figure  2),  which  decreases  the  lateral  spreading 
when  advection  is  present  and  extinguishes  the  signal  when  E  =  0. 

The  effect  of  water  content  on  lateral  spreading  is  illustrated  again  in  Figure  7,  which  is  plotted 
logarithmically. 
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Figure  6:  Effect  of  evaporation  rate  on  the  flux  from  the  surface  as  a  function  of  lateral  position. 
Same  conditions  as  in  Fig.  4  except  that  6  =  0.25.  The  flux  reaching  the  surface  when  E  =  0  is 
negligible. 


Figure  7:  Effect  of  water  content  on  lateral  spreading  of  the  surface  flux. 


Plane  Source  Emission  From  Constant  Concentration  at  the  Depth  of  the  Land  Mine: 
In  this  case 


g{x)  =  Co  0<x<W 
g{x)  =  0  W  <  X  <  L 


Applying  this  to  (35)- (36),  we  obtain 


Eo 


Ej\’  = 


2Co 


N7rZl,r{-P) 


CqW 

'  LZoi-P) 

s.n(— ) 


AT  =  1,2,... 


(43) 


(44) 

(45) 


8 


and  thus  that 


Ct{x,z) 


CoWZojz) 

LZoi-P) 


+ 


2Co 

TT 


Soo 
N=\ 


.Nttx.  ,  ^NttW. 
cos[-^]  sin[— 


NZr,{-P) 


(46) 


J(x,0)) 


CoWQo(O) 

LZoi-P) 


OO 

N=1 


fNirx.  .  .NirW. 

cos[— ^]sm[-Y-] 


Qa^O) 

NZr,{-P) 


(47) 


Figures  8-10  show  compai'isons  of  the  surface  flux  emission  from  a  buried  land  mine  emitting  a 
constant  flux  of  2^  =  10”"^®  g/cm^/ sec,  and  a  zone  of  constant  concentration  Cq  =  4.6  x  10  3  /xg/cm^ 
at  the  same  location. 


Figure  8:  Comparison  of  constant  flux  F  =  10  g/cm^/sec  and  constant  concentration  Cq  = 
4.6  X  10“3  /xg/cm^  sources  at  2  =  —10  cm  for  the  case  of  infinite  half  life  and  E=  0.5  cm/d. 


Figure  9:  Comparison  of  constant  flux  F  =  10“^®  g/cm^/sec  and  constant  concentration  Cq  = 
4.6  X  10~3  /xg/cin^  som'ces  at  2  =  —10  cm  for  the  case  of  r  =  1  yr  half  life  and  E=  0.5  cm/d. 
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Figure  10;  Comparison  of  constant  flux  F’  =  10  g/cm^/sec  and  constant  concentration  Co  = 
4.6  X  10“3  /ig/cm^  sources  at  z  =  —10  cm  for  the  caise  of  r  =  1  yr  half  life  and  E=  0.05  cm/d. 


Continuous  Point  Source  of  Flux  Emission: 

As  the  width  W  of  the  land  mine  is  decreased  toward  zero,  while  the  product  E  =  FW  is  held 
constant,  we  approach  the  delta  function  solution 

f(x)  =  ES(x)  (48) 


which  results  in 


^  2Zo(2)  ,  NTTx  Zjs^iz) 

J(x,0))=|g^  +  2EE^,.iCos[^’^^^^^ 


(49) 


(50) 


LQo  - - 

Figure  11  shows  a  comparison  between  the  point  source  emission  solution  and  one  for  which 
W  —  15  cm  and  F  =  10“^®  g/cm"/sec. 


Figure  11:  Comparison  of  delta  function  (W  «  0  cm)  and  distributed  (W=15  cm)  constant  flux 
sources  at  z  =  —10  cm  for  the  case  of  r  =  1  yr  half  life  ,  E=  0.5  cm/d,  and  E  =  FW  = 
1.3  X  10"'^  /xg/cm/d. 
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